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XXVIII. The Propagation of Heat in Liquid Helium 
at Temperatures below 0-6°K 


By D. V. OsBornet 
Royal Society Mond Laboratory, University of Cambridget 


[Received August 10, 1955 ; revised November 15, 1955] 


ABSTRACT 


Experiments are described in which sine waves of heat are propagated 
in liquid helium at temperatures between 0-25°K and 0-6°x. Although 
propagation takes place in a closed tube, resonance does not occur, but, 
by measuring the phase and amplitude of the temperature waves, values 
are deduced for the velocity and attenuation of the waves at frequencies 
between 60 c/s and 2ke/s. The velocity increases with frequency ; at 
0-5°K it tends to a limit of about 133 msec! at high frequencies, but 
at 0-25°K it reaches about 190 m sec! at 2 ke/s and is still increasing. 
The attenuation per wavelength decreases with increasing frequency 
being nearly 27 at the lowest frequencies. 

At temperatures up to 0-4°K it is shown that the mean free path of 
the phonons must be greater than the tube diameter (2cm) and that 
their reflection at the wall must be partly diffuse. The data are 
quantitatively explained if about one-third of the phonons suffer diffuse 
reflection and the remainder suffer specular reflection. 

Above 0:5°K the mean free path is found to be smaller than the tube 
diameter, being about 0-08 cm at 0-54°K and about 0-04 cm at 0-58°K. 


§ 1. INTRODUCTION 


Stnce the original prediction of Landau (1941) that the velocity of 
second sound in liquid helium II would rise to 1/\/3 of the velocity of 
ordinary sound as the temperature approached absolute zero, several 
experimental investigations have been carried out to test his theory 
(Peshkov 1948b, Maurer and Herlin 1949, Pellam and Scott 1949, 
Atkins and Osborne 1950, Mayper and Herlin 1953, Kramers, van den 
Burg and Gorter 1953, de Klerk, Hudson and Pellam 1953, 1954, 
Kramers, van Peski-Tinbergen, Wiebes, van den Burg and Gorter 1954). 
In all this work (except that of Peshkov whose measurements extended 
down only to 1:03°K) the pulse method has been used and two main 
results have emerged: that at temperatures below 0-6°K the received 
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pulse is greatly rounded and lengthened compared with the input pulse, 
and that the time of arrival of the beginning of the received pulse 
corresponds to a velocity of propagation of 130 to 200m sec’. Since 
the probable velocity of ordinary sound at absolute zero is 239 m sec™? 
(Chase 1953, Chase and Herlin 1955), while the velocity of second sound 
at 1°x is only about 20msec~4, Landau’s prediction is qualitatively 
verified and with it the hypothesis that phonons are the dominant 
excitations in liquid helium at temperatures below 0-6°K. It remains 
to explain the elongation of the received pulse and to explain why the 
velocity of propagation tends to exceed the Landau value of 138 msec™? 
at the lowest temperatures. Discussions of these phenomena in terms 
of the behaviour of the liquid have been given by Dingle (1952), Osborne 
(1952), Gorter (1952), Atkins (1953), Mayper and Herlin (1953), and 
Kramers et al. (1954), the last of these papers containing also some 
discussion of the behaviour of the heater and thermometer. 

In the present experiment sine waves are used so that the pulse 
distortion previously observed appears as measurable dispersion and 
attenuation. Moreover the method is such that the velocity and 
attenuation measurements do not involve any assumption concerning 
the behaviour of the heater and thermometer at high frequencies. 

The sine waves are emitted by a plane heater at one end of a closed 
cylindrical tube and detected by a plane thermometer at the other end. 
This technique was first used by Peshkovy (1948 a), who observed 
resonance in the tube and so deduced the velocity of the waves at 
temperatures above 1°K. It is now found that no resonance can be 
observed at temperatures below 0-7°K, but that it is still possible to 
estimate the velocity and attenuation of the waves by measuring the 
phase and amplitude of the sine waves received at the thermometer. 
The use of sine waves means that a detector of very narrow bandwidth 
and correspondingly low noise level can be used, and this facilitates the 
measurement of the very small signals obtained in the absence of resonance. 

Since the thermometer has an unknown frequency response and phase 
lag, experiments are carried out with several different lengths of tube 
using the same heater and thermometer each time. A final run, using 
the first tube length again, is necessary to verify that the thermometer 
behaviour does not change with time. The velocity and attenuation at 
any given frequency are then obtained by comparing observations made 
at that frequency during each of the five runs. The experiments are 
described in §§ 2 to 5 and their interpretation in §§ 6 to 9. 


§ 2. Tur Cryostat 


The problem of cooling liquid helium by adiabatic demagnetization is 
complicated by the two opposing requirements that it must be easy to 
dispose of the heat of magnetization, but that the heat leak to the 
helium, once cold, must be as small as possible. The device used for 
solving this problem in the present cryostat is due to Dr. J, Ashmead 
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and one use of it has already been reported by Herlin and Chase (1955) : 
it is illustrated in fig. 1. The experimental chamber A is surrounded ne 
a vacuum space and is connected to the helium bath B by a conical 
German silver tube C of length 6 cm, wall thickness 0-3 mm, diameter 
at the larger end 1-6 cm and diameter at the smaller end viva Into 
this tube fits a plug D consisting of a similar cone closed at both ends 


Fig. 1 
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The cryostat. L, helium dewar and liquid air shield; for remaining key see § 2. 


and evacuated. The plug may be raised or lowered by a control rod 
projecting out of the top of the cryostat. The experimental chamber 
contains the paramagnetic salt E, which acts as the refrigerant, and 
the experimental apparatus F, both immersed in liquid helium. Initially 
the plug is raised so that there is direct communication between the 
experimental chamber and the main helium bath at 1-2°K and with the 
apparatus in this condition a magnetic field is applied to the salt by 
means of the magnet G, an armoured water-cooled solenoid with its 
axis vertical. The heat of magnetization is led away at once by the 
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helium in the conical tube to the helium bath and, on account of the 
excellent thermal conductivity of liquid helium at 1-:2°K, temperature 
equilibrium is reached in a few seconds. The plug is now lowered so as 
to fit tightly in the conical tube, and the magnetic field is reduced to 
zero. The experimental chamber cools, in the present instance to 0-2°K, 
and may take between 20 and 40 minutes to return as far as 0-7°K, the 
warming rate depending partly upon the amount of power supplied 
during the actual measurements. The conical plug and tube have 
well-machined surfaces and, although it is easily possible for superflow 
to take place past the plug, the fit is sufficiently good to offer considerable 
resistance to the passage of the normal component of liquid helium, 
the flow of which is hindered by an ordinary viscosity of the order of 
20upoise. Since it is the normal component which carries heat, the 
heat leak is by this means reduced to quite a small value. 

The use of a movable plunger in this way has the advantage over the 
more usual capillary tube that large heats of magnetization can be quickly 
conducted away. It is therefore particularly suitable for experiments in 
which a large mass of refrigerant salt has to be used, a condition which 
obtains in this experiment both because of the large volume of helium to 
be cooled and because of the occasional necessity of using relatively large 
heat inputs for experimental purposes. The refrigerant E in the present 
work consists of 200g of coarsely powdered manganous ammonium 
sulphate placed above the experimental space so that it may absorb 
the residual heat input from the cone, leaving the region below relatively 
free from unwanted heat currents. 

The experimental space has a capacity of 150 cm’. Electrical leads 
for the heater and thermometer are introduced via two 1mm tubes 
packed with petroleum jelly, one of the tubes being shown at H. These 
almost certainly allow the passage of superfluid from the experimental 
chamber to the helium bath but, like the cone, they do not permit a 
sufficient flow of normal fluid in the reverse direction to constitute an 
appreciable heat leak. 

In an appendix below the experimental space is a small ellipsoid J of 
copper potassium sulphate which is the core of a mutual inductance K 
and thereby enables the temperature in the chamber to be measured. 
This ellipsoid is never magnetized and is cooled solely by thermal contact 
with the helium surrounding it, so that it cannot attain a temperature 
lower than that of the helium. In fact it always responds within two 
or three seconds to any change of temperature produced by changing the 
power dissipated in the second sound apparatus. Since the rate of warming 
is not more than 0-01 deg. min-?, this means that the thermometer gives 
a faithful indication of the helium temperature. The salt J is that 
originally used by Garrett (1950), and the absolute temperature is deduced 
as indicated in his paper, due allowance being made for the magnetic 
behaviour of the German silver tubes. The necessary correction was 
estimated in some auxiliary experiments. 
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About 6 litres of liquid helium are required to cool the cryostat from 
liquid air temperatures down to 4-2°K and to fill it. The length of run 
is then about five hours of which one hour is taken up in calibrating the 
measuring salt J at liquid helium temperatures. During the remaining 
time, four or five demagnetizations are carried out, the actual process of 
demagnetization taking about five minutes and the warming time between 
0-2°K and 0-7°K being about forty minutes. Ina trial run it was established 
that the natural heat leak was of the order of 600 erg sec-! and that 
the rate of rise of temperature varied from about 2 mdeg. min-! at 
0-25°K to about 7 mdeg. min~ at 0-5°k. The thermal capacity which 
may be deduced from these values is in reasonable accord with that 
estimated from the known masses of salt and helium. 


$3. THe PROPAGATION TUBE 


That part of the second sound apparatus which is in the cryostat is 
shown in fig. 2. The tube A defining the separation between the heater 
and thermometer is of silica and has an inside diameter of 2:02 cm and 
a wall thickness of 2mm. The heater B consists of some 70cm of 


Fig. 2 
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46 s.w.g. Eureka wire wound zig-zag on a thin paxolin former and 
backed by a glass disc C 3mm thick. The thermometer D is a thin 
strip of mica, 2cmx2mm x20, painted over with Aquadag and 
silvered at the two ends for better electrical contact. It is laid across 
a diameter of the tube with its coated side facing the heater, and the 
silvered ends press against copper contacts ER’ Sule to the silica with 
Sellotape. It is backed by a second glass disc F. The whole was held 
together with Sellotape during assembly and was finally clamped 
together by a light copper frame. There are sufficient slight gaps at the 
ends to admit helium freely and to maintain temperature equilibrium with 
the rest of the inner chamber. Current and potential leads for the 
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heater and thermometer are led up the two tubes H (fig. 1) to the 
helium bath, the heater leads being kept separate from the thermometer 
leads to reduce the coupling between them when using sine waves. 

When the heater is switched on to generate second sound, the 
temperature in the tube rises with a time constant of about two seconds 
to a value slightly higher than the temperature measured by the salt. 
The effect is negligible except at the lowest temperatures when using the 
highest powers. In these circumstances an appropriate correction has 
been applied to give the temperature inside the chamber from the 
measured temperature of the salt J. The final accuracy of the temperature 
assigned to any given phase and amplitude measurement is -—-10~? deg. 


§ 4. THe ELECTRICAL EQUIPMENT 


Outside the cryostat are the oscillator which supplies sine waves to 
the heater, and the amplifier and cathode-ray tube which amplify and 
display the signals from the thermometer. The amplifier incorporates a 
phase sensitive detector which makes the bandwidth of the system about 
0:2c/s and the resulting noise level is sufficiently low for input signals 
of the order of 0-1 xv to be detected. The cathode-ray tube is so connected 
that the polar coordinates (v, «) of its spot give a direct measure of the 
amplitude and phase angle of the amplifier output. 

In order to calculate the amplitude and phase of the thermometer 
waveform from the observed values of v and e¢, the gain and phase shift 
of the amplifier system have been measured at all relevant frequencies. 
The temperature amplitude may then be calculated by dividing the 
amplitude of the fluctuation of thermometer resistance by the slope of 
the resistance-temperature calibration curve but, except at the lowest 
frequencies, this procedure is unlikely to give the correct absolute 
values because the steady state calibration is unlikely to be applicable 
at audio frequencies. Similarly the phase of the observed fluctuations 
of thermometer resistance is likely to lag behind the phase of the 
temperature fluctuations in the helium except at very low frequencies. 
Since, however, the same thermometer was always used, it is legitimate 
to compare the effective amplitudes and phases obtained in this way at 
one given frequency and temperature for different lengths of propagation 
tube. 

The accuracy of phase measurement is better than +2 degrees for 
sufficiently large signals but for very small signals (particularly at high 
frequencies, at higher temperatures and with long tubes) the accuracy 
falls in some cases to +-5 degrees. Fortunately it is in these circumstances 
that the phase shift is large and the accuracy less important. 


§ 5. EXPERIMENTAL RESULTS 


The whole experiment consisted in carrying oui four successful runs, 
each with a different tube length but with the same heater and thermo- 
meter, and then repeating the first run to check the stability of the 
thermometer. The tube lengths used were 0-38 em, 3:16 cm, 9-26 em 
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and 12:35cm. In order that the readings of different runs might be 
compared, observations in each run were made at fifteen fixed frequencies 
(f) between 60 c/s and 2ke/s. For each frequency the polar coordinates 
(v, <) of the cathode-ray tube spot were plotted against temperature 
and the readings corresponding to six fixed temperatures between 0-25°K 
and 0-58°K were read from each curve. Readings at a given frequency 
and temperature could then be compared for different runs, i.e. for 
different tube lengths. A representative selection of the readings is 
shown in figs. 3 and 4. 

In fig. 3, where phase, ¢, is given as a function of tube length, L, for 
the various frequencies, ¢ has been obtained by subtracting the known 
phase shift of the amplifier from the observed phase, «. There remains 
the unknown phase lag of the thermometer, a correction for which would 
decrease the values of 4 particularly at the higher frequencies (see § 8). 
Points in fig. 3 corresponding to a given frequency have been joined by 
straight lines merely to facilitate identification. { 

The observed amplitude, v, is proportional to 0, the temperature 
amplitude of the received signal and @ is calculated from v using the 
known gain of the amplifier and the known sensitivity of the thermometer. 
That this sensitivity is strictly applicable only to steady temperatures 
means that the values of @ so calculated will in general be smaller than. 
the true temperature amplitude in the helium, especially at the higher 
frequencies (see § 8). The values of @ are divided by the amplitude, W, 
of the power flow per unit area and it is found that for the given values of 
frequency and tube length, all the values of 6/W fall on a smooth curve 
if plotted against temperature, ie. that under fixed conditions of 
temperature, frequency and tube length, 9 must be proportional to W. 
In presenting the amplitude data, therefore, it has been found convenient 
to give the quantity A, defined by 


A=6/Wp. 


y) is here the characteristic impedance of liquid helium for second sound 
waves (Osborne 1948) and it serves both to relate the values of A obtained 
at different temperatures and also to render A dimensionless. y) is equal 
to (pc,C)-! where pC is the specific heat per unit volume and c, is the 
velocity of second sound. Since the latter quantity is under discussion, 
the value of y, actually used has been that based strictly on Landau’s 
theory, using the single experimental datum that the extrapolated 
velocity ¢ of ordinary sound in helium at the absolute zero is 239 m sec"! 
(Chase 1953). From Landau’s theory 
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The experimental value of C (Kramers et al. 1952) is about 10% higher 
than the Landau value, and the experimental values of c, so far published 
range from about ¢/1/3 to c. 

The first and fifth runs, both with L=0-38 em, give phases which agree 
to within —-3 degrees. The amplitudes obtained in the fifth run are 
systematically about 15°, less than those of the first run but in view of 
the very wide range of amplitudes used this error is not regarded as very 
serious. The phase shown in fig. 3 for this tube is the mean of the two 
runs ; the amplitude shown in fig. 4 is that obtained during the first run. 


Table 1. Typical Signal Sizes 


Length Frequency | Temperature Power tae 

Lem fe/s Kk W pW em™ 6 udeg. 
0-38 60 0-25 1-0 110 
0-38 60: 0:5 2-9 53 
0:38 2000 0:25 22 oo 
0:38 2000 0-5 22 1-4 

12-35 60 0-25 5:3 19 

12-35 2000 0-5 49 0-27 


Table 1 gives some idea of the orders of magnitude used. For a given 
power the received signal decreased with increasing length, with increasing 
frequency and with increasing temperature. The signal amplitudes 
quoted are those recorded by the thermometer. At 60 c/s these will be 
about the same as the amplitudes in the helium but at 2000c/s the 
amplitude in the helium may be several times larger than that recorded 
by the thermometer (see § 8). 

At temperatures between 0-2°K, the lowest reached, and 0-3°K the 
measured phases depend alittle on the power used. At the lower frequencies 
an increase of power by a factor of ten can produce a phase advance up to 
15 degrees. The phenomenon suggests that second order effects of two 
kinds may become important ; the amplitude of the normal fluid velocity 
may become comparable with the wave velocity (Osborne 1951, Temperley 
1951, Khalatnikov 1952), and the average normal fluid velocity due to 
the mean flow of heat along the tube may become comparable with the 
wave velocity. It may be shown, however, that both of these effects 
might be expected to change the velocity of second sound by a fraction 
of the order of magnitude of 36/7’, and reference to table 1 shows that 
this fraction is in all circumstances very small compared with unity. 
The observations, therefore, remain at present unexplained. 


§ 6. Tue Continuum THEORY 
If the phonons have very short mean free paths, the fluid may be treated 
as a continuum, and the observed amplitudes and phases will be explicable 
in terms of the propagation of waves. having a definite velocity c’, 
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attenuation «, and characteristic impedance y. The observations will 
be treated on the basis of this assumption and it will be shown that the 
result is satisfactory only at temperatures above 0-5°K. It will therefore 
be inferred that below 0-5°K the mean free paths are no longer short, and 
an alternative interpretation of the low temperature data will be given 
in § 7. 

If the characteristic impedance y of the helium for a travelling wave 
propagated in the tube is defined by 
“dg temperature amplitude 

heat flow amplitude per unit area 


and if the thermal penetration impedance ¢ of the two end walls is 
given by 


id 


ne temperature amplitude at end wall 
Y heat flow amplitude into end wall per unit area’ 


then the usual methods of transmission line theory may be used to show 
that the temperature amplitude @ at the thermometer will be given by 
g ¢ 
Wee oh plac -cuaaith pies. ° y - ¥) 
where W=heat flow amplitude per unit area given out by heater, 
L=length of tube, B=a+ 7k, «=attenuation constant per unit length, 
k=2n/A and. A=wavelength. y and ¢ will in general be complex. 

If the reflection coefficient is very nearly unity (| ¢|>1) we shall 
expect @ to be inversely proportional to L, at least for the two shorter 
tubes, and this is in fact observed at the lower frequencies, as may be 
verified by reference to fig. 4. 

The very high value of | ¢| arises partly because the end plates are 
of thermally insulating: material, and partly because of the Kapitza 
(1941) boundary resistance which impedes the flow of heat from the 
helium into the solid. 

At the higher frequencies, where even the 3:16 cm tube becomes a 
considerable fraction of a wavelength, L will no longer be sufficiently 
small for this approximation to hold, and it is indeed observed (e.g. fig. 4(0), 
curve F) that the inverse proportionality breaks down. At these 
frequencies no simple argument is possible to show that | ¢|>1, but 
the assumption that this is so is justified on the one hand by the fact 
that ¢ would not be expected to vary very rapidly with frequency and on 
the other hand by the fact that the interpretation of the high frequency 
results is relatively insensitive to the value of ¢. 

Upon making the assumption |¢ |>1, we find from eqn. (1) a 
theoretical amplitude A’ and phase lag ¢’ given by 


‘ @ , ny  & \) -1/2 5 

A'= ow = {sin +-sinh a) fi Re ak 5) 

and tan ¢’=tan y coth G #) oa CP ORCS SLE) 
where p=kL=2rfL/c’, 
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o'(s) and A’(xs) for a few representative values of «/k are shown in 
figs. 5 and 6 respectively. For convenience in the use of tables, «/k 


has been expressed in radian deg.~?. 


Fig. 5 
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Theoretical phase lag. Valuesof «/k in radian deg.—! are indicated for each curve. 


The experimental points ¢(Z) and A(L) (figs. 3 and 4) are fitted to 
these curves by plotting 4(L), A(L), ¢/(%) and A’(y%s) on transparent 
logarithmic graph paper and superposing each set of experimental points 
in turn on the appropriate theoretical curve. «/k is determined by 
selecting that theoretical curve which fits the experimental points best, 
and the relative position of the abscissae ;s and L is adjusted for the best 
fit. This adjustment determines the constant of proportionality relating 
ys to L and hence gives the velocity c’. In order to allow for the unknown 
phase lag and efficiency of the thermometer, the relative position of the 
ordinates is also adjusted for the best fit. At any given temperature and 
frequency the same pair of values of «/k and c’ must, of course, be used 
to fit both the phase and the amplitude data. It proves possible to carry 
out this fitting process within the limits of experimental error for all 
frequencies and for all temperatures, except for one or two individual 
points which seem rather far out, 
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The results are given in figs. 7 and 8 where the wave velocity c’ and 
the attenuation per wavelength «A (=360a/k if a/k is in radian deg.) 
are given as functions of frequency at various temperatures. ‘The accuracy 
of the velocity determinations improves with increasing frequency, the 
values being accurate to within about +20% at 60 c/s, +5% at 480 c/s 
and +3% at 2kc/s. The attenuations are probably accurate to within 


Fig. 6 
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Theoretical normalized amplitude. 


Values of «/k in radian deg.—! are indicated for each curve. 


+20°%, at all frequencies. The chief source of error at low frequencies 

is the possibility of fitting the data over rather a large range of values 

of «/k and the chief source of error at high frequencies is the smallness 
ignals. 

re eater these velocities and attenuations by treating the 

normal component of liquid helium II as a continuous medium with a 

definite shear viscosity, second viscosity, and thermal conductivity. 
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We have then the following equations of motion for the normal component 
and for the superfluid : 


wee Vi es 5 (VNataV(V v,))— on V(V -¥q)—=0 . (4) 
nD 


and AL 8 + Vp—SVP=0, (5) 


where v,, and v, are the velocities of the two components, pn and p, are 
their densities, is the viscosity of the normal component as measured 
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Wave velocity c’ as a function of frequency. 


(a) T=0-25°x, (b) T=0-3°x, 


(c) T=0-4°x, 
(d) T'=0-5°x, 


(c) T=0-54°K, (f) T'=0-58°x, 
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(d) T=0-5°K, (ec) T=0-54°x,  (f) T=0-58°K. 


for example, in a rotating cylinder viscometer, jw, is the second viscosity 
of the normal component and S is the entropy per unit mass of helium. 
The equations are valid only for v, and v, small compared with the velocity 
of sound, and for small excursions of p, p and 7’ from their mean values. 
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The equations of continuity of matter and continuity of entropy are 


a 
ay EV - (Pa¥at Pave) —=0 | ge 

A(pS) K 
and tee) + pSV. va— 7 VEN=0 _'. 23 ees 


to the same approximation, where K is the thermal conductivity of the 
normal fluid. We shall assume that the thermal expansion coefficient 
(Atkins and Edwards 1954, 1955) is very small so that p may be taken as 
a function of p only, and S as a function of 7’ only. For second sound 
above 1°K where dissipative processes are usually unimportant, the net 
flow of matter is zero, i.e. py¥,--ps¥;=9. This is not necessarily true 
when dissipative forces are acting, but we shall assume that ppp and 
p,v, are of the same order of magnitude and accordingly, since p,> py 
we have v,<v,. This is equivalent to considering the behaviour of a 
phonon gas in a fixed medium (see Ward and Wilks 1952). Equations (4) 
and (7) may then be treated by the method of Kirchhoff (1868) (see also 
Rayleigh 1896, Weston 1953). (For an alternative treatment following 
Helmholtz 1863, see Dingle 1950.) 


Table 2. The Continuum Theory 


x < ‘Narrow tube’ x > ‘ Wide tube’ 
¢ 4/x x tf 1 1 19 
= 3 ( gt ani) (1 = s/n) eae) aaa) 
2 2 1 ul 

d Sat Bb =) a AT® [ae ee oe 
n( 6° 73 V (aa) \* a/Qa) + 4(2m) 

7 12; he nin 
ny, ore {ae ant 
Yo f/x % V (22) (1+ Bez) 


Two approximate solutions emerge, depending upon the ratio of 8 to a, 
where 6 is the penetration depth of viscous waves given by 8?=p/zp,f 
and a is the radius of the tube. If 5/a<1 only a thin layer of the medium 
near the wall is affected by viscosity, and the attenuation and change of 
velocity are small; but if d/a>1 the flow reduces to an alternating 
Poiseuille flow, and the propagation characteristics are profoundly 
modified. The condition for the validity of the calculations is that the 
mean free path / shall be short compared with each of the four lengths §, 
a, X and L, the last two being the wavelength and the length of the tube 
respectively. It is assumed, in accordance with the usual results of kinetic 
theory, that the thermal diffusivity (K/pC), kinematic viscosity (/p,) and 
kinematic second viscosity (jo/p,) are all of the order of magnitude cl 
where ¢ is the velocity of phonons. With these approximations and 
assumptions, thermal conduction and second viscosity disappear from the 
final expressions, 


Cleo 


0-5 


0:3 


0:2 
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The results for velocity and attenuation, given in table 2, are expressed 
in terms of the dimensionless parameter x=27fa2/u’, where p’ is the 
kinematic viscosity, /p,. > 1 corresponds to §<a and x <1 corresponds 
to >a. The velocity c’ is given in terms of the velocity c, of second 
sound in the free liquid in the absence of dissipative effects. 7 is the 
effective value of thermal impedance defined by y=0/W. W is the heat 
flow amplitude averaged across the whole area of the tube on the 
assumption that the heater is a uniform plane source. @ is the temperature 
amplitude and this is averaged along a diameter of the tube to correspond 
with the fact that the thermometer is a thin strip across one end. yg is 
defined as (pC'c,)~1 where pC is the specific heat per unit volume, and is 
therefore the thermal impedance of the free liquid in the absence of 
dissipation. The expressions for c’/c, and for «A contain higher order 
terms in x and these have been obtained by following Weston’s (1953) 
extension of Kirchhoff’s treatment for ordinary sound. 
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Theoretical curves of velocity against x. The two parts of this curve 

represent the approximate formulae of table 2 for small x and for 


large x. 


Logarithmic graphs of the functions c'/c, and «A are given in figs. 9 
and 10, Since « is proportional to frequency these curves may be compared 


Y2 
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with the logarithmic curves of c’ and «A against f (figs. 7 and 8) it and will 
be seen that they are of the same general form. In order to make a detailed 
comparison the method of superposing logarithmic graphs is again used, 
the values of «A as a function of frequency (fig. 8) being first superposed 
on the theoretical curve of «A versus a (fig. 10). The theoretical and 
experimental values of the ordinate «A have to agree without any adjustable 
constant of proportionality, and the theory must therefore be made to fit 
by adjusting the relation between the abscissae. This process determines 
the ratio of x to f, ie. it determines the quantity 27a?/u’ and hence 


Fig. 10 
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Theoretical curves of attenuation per wavelength against x. The two parts 
of this curve represent the approximate formulae of table 2 for small a 
and for large x. 


determines y’. It is found that a tolerable fit can be obtained at all six 
temperatures by the choice of an appropriate value of f/x at each 
temperature ; the high frequency results are used to determine f/x and 
the low frequency results are, within rather large experimental errors 
also consistent with the theory. The values of «’ obtained in this war 
(using the fact that a=1-01 cm) are shown in table 3. In order to 
determine c,, the experimental curve of c’ versus f is superposed on the 
theoretical curve of c'/cy versus x, making use of the relation between f 
and w already found. When the ordinates are suitably adjusted it is found 
that a good fit can be obtained, again within rather large experimental 
errors at the low frequency end. The value of c’ on the experimental 
graph which coincides with the ordinate c’/c,=1 on the theoretical graph 
is taken as the value of c, and is given in table 3, ; 
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The actual viscosity j4 may be obtained from the kinematic viscosity ju’ 
if p, is known. p, may be obtained from the Landau theory using the 


extrapolated velocity c of ordinary sound for the velocity of phonons ; 
this procedure gives 


Pu= 3p l'C/e?=1-78 x 10-574 g cm-3. 


We may also make use of the relation derived from elementary kinetic 
theory, 


p= sel, 
to calculate formally a kind of mean free path J for the phonons. Values 
of p,, » and J are given in table 3. If the continuum theory is to be 
self-consistent, this mean free path must be shorter than each of the 
four lengths 6, a, A and L. We note from table 3 that / is always less than 
these lengths, and in most cases less by a considerable factor. 


Table 3. Results according to the Continuum Theory 


Kaneratie Normal Mean _ | Penetration 
Temp ; ; Velocity fluid Viscosity | free depth at 
oes viscosity wa : ; : 
ORM, _3| Co msec density ppoise | path 2 ke/s 
po cm*sec—} 3 ie 
Pag cmt tom 5 cm 
0-25 1025 210 6-95 x 10-8 val 0-13 0-40 
0:3 830 180 14-4 10-8 119 0-11 0-36 
0-4 740 160 44-5 10-8 320 0-092 0-34 
0:5 545 140 Ti. 10-% 600 0-068 0-30 
0-54 640 145 15.10% 970 0-080 0-32 
0-58 320 115: 202 x 10-8 650 0-040 0-23 


The continuum theory, then, appears to offer a self-consistent description 
of the experimental facts at all temperatures from 0-25°K to 0-58°K. It 
has however two major shortcomings at the lower temperatures: it 
demands the assumption that c, increases to as much as 210 msec 
at 0-25°K, and it gives a viscosity which increases with increasing 
temperature between 0-25°K and 0-54°K. These two results are in 
serious conflict with the theoretic expectations of Landau (1941) and 
of Landau and Khalatnikov (1949) respectively. It is perhaps particularly 
noticeable that the mean free path / given in table 3 diminishes by less 
than a factor of 2 when the temperature rises from 0:25°K to 0-5°K, a 
rather improbable situation. 

It is therefore believed that the continuum theory may be legitimately 
applied only at the two highest temperatures used here, Le. 0-54°K and 
0-58°k. The conclusions summarized in table 3 are believed to be valid 
only for these two temperatures, and the figures shown for lower 
temperatures are without real physical significance, except that the mean 
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free paths at lower temperatures are probably related to the tube diameter 
and to the degree of specular reflection of phonons at the walls. This is 
discussed in § 7. 


§ 7. Tue Lona Mpan FREE Pato THEORY 


It remains to show that the observations at temperatures below 0-5°K 
may be explained by assuming the mean free path of the phonons to be long 
compared with the tube diameter. Ziman (1954) has treated the problem 
of the propagation of a pulse of heat along a cylindrical tube, assuming 
a mean free path long compared with the diameter of the tube and diffuse 
reflection of phonons at the walls. His treatment gives a correct 
description of the general shape of the received pulse but it cannot claim 
to give a valid description of the beginning of the pulse because it depends 
on the assumption that every phonon makes many collisions with the 
wall during its journey from heater to thermometer. This leads to a 
diffusion equation which in turn predicts that the thermometer temperature 
will start to change immediately the heat pulse is applied. In practice, 
the earliest phonons to arrive at the thermometer will do so after a time 
of the order of L/c by taking comparatively direct paths and will not fall 
within the scope of the above assumption. Moreover, in terms of sine 
waves, this theory is not applicable at high frequencies where the periodic 
time is short compared with the time taken for a phonon to travel from 
one collision with the wall to another. In the absence of a treatment valid 
for all frequencies, however, we shall in the first place adapt Ziman’s 
treatment for use with sine waves and apply it to the low frequency 
results. 

In the steady state the heat flow F along a cylinder whose radius is 
small compared with the phonon mean free path is 


dT 
F=¥ra'pCc— 3 2 ee ee 


(Casimir 1938, with the correction pointed out by Berman, Simon and 
Ziman 1953), implying an effective thermal conductivity which may be 


written in the form p 
Kp pCcl(2a). aos, 1 at oe eee 


where the diameter 2a acts as the mean free path. 

For subsequent comparison with the continuum theory we note that 
a flow of heat in liquid helium II corresponds to a flow of normal fluid 
and that a temperature gradient corresponds, by virtue of the fountain 
effect, to a pressure gradient. Using the standard formulae of the two 
fluid theory, Landau’s formula for the phonon contribution to p,, and 
the fact that the specific heat is proportional to 73, eqn. (8) can be 
transformed to 
27a? dp 


PR Pra? Same Vara cr SOY Gen, 


C= 
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where Q is the volume flow of normal fluid. This relation may be 
compared with Poiseuille’s formula for flow through a tube and it suggests 
an effective viscosity 


1 
Pica ay Puce ag tO PAR he. ep (L1) 


where the effective mean free path is now 3/32(2a). 

At this point we introduce a coefficient €, being that fraction of the 
phonons which are scattered diffusely at the walls, the remaining 1—é 
suffering specular reflection. We also introduce a factor to allow for the 
fact that the tube is not infinitely long and that there is some probability 
of a phonon travelling straight from one end of the tube to the other 
without colliding with the wall. For steady flow these two factors appear 
in eqn. (8) in the form 


77 2 3h Er 
Pa }reepCess W (=), a= Mead pert 12) 


where L'=€L/(2—€). W has been the subject of analysis by Clausing 
(1928) and it is discussed by Kennard (1938) ; it takes the value (Clausing 
1928, p. 103, col. 4) 8a/3L’ for very long tubes. L’ is used here in place 
of L because in Clausing’s work the scattering at the wall has been taken 
to be perfectly diffuse, whereas a proportion of specular reflection reduces 
the effective length of the tube. Using eqn. (12) for # we obtain an 
effective thermal conductivity of 


Kg=4pCet2a) 52 W (=) Pe te ee ah) 


8a 
and an effective viscosity of 
1 8a 1 


Hett= 35 Pal(24) oF WLTa)’ (14) 

The effective mean free path for thermal conduction becomes 

3L Lt 
oS se ==} 7 BOL Seder chi et LD 
Late = w (=) 20 ( 5) 
and that for viscosity becomes 
| 8a 1 3 

Ee ew me (MI gg ve (LG 
le 3, W(Lia) 32°") ue 


Provided, now, that the flow in the tube is substantially uniform. 
ie. provided that the tube is considerably shorter than a wavelength, 
we may substitute the modified thermal conductivity (eqn. (13)) in the 
classical equation of heat conduction, 

athe Vai b 
B= 0 OF 
and obtain the velocity and attenuation of temperature waves as 


a VL, Tip 1/2 
GO = a eee w(=) 


4ca 3L LEN\ v2 
and RN AT Ba W Ge ; 


X 
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Since this is a wave characterized by «A=27, the phases and amplitudes 
predicted by it would be those given by the curve «/k=0-0175 radian deg.~* 
in figs. 5 and 6, were it not for the fact that the velocity and attenuation 
vary with the tube length. This fact necessitates the separate calculation 
of the predicted phase shift and amplitude as a function of tube length 
for each frequency used and for a range of possible values of €. These 
curves are then compared with the experimental points (figs. 3 and 4) 
and the value of € is chosen which gives the best agreement between 
experiment and theory. In making the comparison, the relation between 
the abscissae is fixed because both the theoretical and experimental points 
are plotted against actual tube length in centimetres ; the ordinates may 
however be adjusted to allow for the unknown phase lag and efficiency 
of the thermometer. For frequencies from 60 c/s to 361 c/s, the experi- 
mental amplitudes agree with the theoretical ones within the limits of 
experimental error, but the errors are such that € cannot be determined 
from these measurements alone. The phase data are however much more 
accurate and the theoretical curves are more sensitive to the value assumed. 
for €. It is possible to obtain an estimate of € at each of the six lowest 
frequencies, and the averaging of these gives the values shown in table 4 
where the errors given are the mean deviation of the mean value of € 


Table 4. Fraction of Phonons Diffusely Reflected 


fags é 

0-25 0-26+0-01 
0-3 0-29 +0-02 
0-4 0-35 +0-03 
0-5 0-28 +0-03 


The mean wavelength of phonons is of the order of 200A at 0-5°K, and 
so may well be comparable with the smoothness of the walls. Values of € 
of the order of 0-3 are therefore not unreasonable, but since the phonon 
wavelength decreases with increasing temperature it would be expected 
that the proportion of diffuse scattering would also increase with 
temperature. For this reason it is believed that this calculation is 
probably not applicable at 0-5°K because the mean free path is not 
sufficiently long, although it is thought to be applicable between 0-25°x 
and 0-4°K, 

For reasons already mentioned the foregoing treatment is valid only for 
low frequencies ; an idea of the limiting behaviour at high frequencies 
may be obtained by the following device, which is probably capable of 
giving accurate predictions if suitably developed. If we consider the 
emission from the heater of a very short pulse of heat, the temperature 
fluctuation @ at a point distant z from the heater would be given, according 
to Kramers et al. (1954) following Ziman (1954), approximately by 


O(7)=Kxr-V2 exp (— 1/7) 
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where T= 8act/32° and « is a constant. This cannot be strictly true since 
it predicts a temperature which starts to rise immediately the pulse is 
emitted, whereas there must in reality be a delay of z/c before any phonons 
can reach the point of observation and before any effect can be observed. 


As a trial function, we shall assume that a pulse would in fact give rise to 
a fluctuation 


O(7)=«(r—79)¥? exp {—1/(r—79)} ; a a 


O(7r)=0 STS.) 


where t9=(8ac/32")(z/c). We then regard sine waves of heat emitted 

from the heater as a succession of pulses of varying amplitude. By adding 

together the effect of all such pulses preceding the instant of observation 

ne obtain for the resultant temperature fluctuation at distance z from the 
eater 


Sac \ 1/2 8w \ 1/2 38w \ V2 @ 
ee) af) 2 g 
(t, 2) (¥) exp Ga) psn {(Z2) abe wt- aah 


(18) 


At low frequencies the term wz/c may be ignored, and the wave velocity 
and exponential attenuation agree with those previously found. The 
factor z-! appears to be spurious and has probably arisen because the 
trial function is not sufficiently accurate, but the expression as a whole 
works sufficiently well at low frequencies to engender some faith in its 
predictions for high frequencies. For frequencies so high that 


w/e> (3w/4ac) V2, 


the wave velocity is predicted to be simply c, and the transition frequency 
at which the two terms have equal importance is /j>=3c/87a=2:8 ke/s. 
fp is such a frequency that its period is comparable with the time taken 
for a phonon to cross the tube. 

The foregoing is to be regarded as strong, although not conclusive, 
evidence that, if the phonons have a long mean free path and suffer 
partially diffuse reflection, the velocity of temperature waves will tend to 
that of ordinary sound at very high frequencies. The curves of fig. 7 
show such a tendency at temperatures of 0-4°K and below, but not at 
higher temperatures. The derivation of these curves is of course only 
strictly applicable if the continuum theory applies because the end effects 
become important if the mean free path is long, and the velocity and 
attenuation then depend on length. It is clear also, however, from direct 
inspection of the phase difference between the two longest tubes (for 
which the end corrections are smallest), that the rate of change of phase 
with length certainly involves a velocity greater than 140 sec"! at 
frequencies and temperatures where the curves of fig. 7 exceed this value. 
From the results in the whole frequency range we conclude, therefore, 
that the mean free path is substantially longer than the tube diameter at 
temperatures of 0-4°K and below, but that at and above 0-5°K this is no 


longer true. 
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It is to be noted from the earlier part of the above analysis that at low 
frequencies the ‘ open gas’ of phonons flows as though it were controlled 
by a viscosity determined by eqn. (14), implying an effective mean free 
path given by eqn. (16). With the values of € given in table 4, this 
effective mean free path is about 0-1 cm for the longer tubes. This rather 
small value accounts for the deceptive self-consistency of the continuum 
theory when applied at low temperatures where the mean free path is in 
fact long (see table 3). The mean free paths of table 3 at lower temperatures 
cannot, however, be quantitatively compared with eqn. (16) because they 
were deduced from the data at all frequencies whereas eqn. (16) is 
applicable only to low frequencies ; it can only be said that the agreement 
in order of magnitude is satisfactory. 


§ 8. THE THERMOMETER RESPONSE 


Throughout the calculations allowance has been made for the fact 
that the phase and amplitude recorded by the thermometer may differ 
from those existing in the helium by amounts which will depend on 
frequency and on temperature but not on tube length. Having established 
the theoretical interpretation without assuming any values for these 
unknown quantities, it is now possible to calculate the actual amplitudes 
and phases which ought to be observed and, by comparing them in 
absolute magnitude with those actually found, to deduce the phase lag 
and the efficiency of the thermometer. The latter is defined as the ratio 
of the temperature amplitude given by the thermometer to that actually 
existing in the helium. 


Table 5. Thermometer Response 


At 60 ¢/s At 361 c/s At 2 ke/s 
shete has: Phase Phase Phase 
lag | Efficiency | lag | Efficiency lag | Efficiency 
(deg.) (deg.) (deg.) 

0-25 28 0-8 88 0-14 

0-3 14 0-9 82 0:3 

0-4 4 0-9 64 0-5 

0-5 0 0-8 123 0-4 

0-54 0 0-8 125 0-4 

0-58 0 (1:1) 118 0-8 


The result of this process is shown in table 5. The long mean free path 
theory which is applicable at the lower temperatures can be applied 
quantitatively only at frequencies up to about 361 c/s. At the higher 
temperatures, the continuum theory can readily be applied at 60 c/s 
and at 2ke/s, but the calculation of absolute amplitude and phases 
presents some difficulty in the middle frequency range where the tube 
is neither ‘ wide’ nor ‘ narrow ’. 
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The figures in table 5 apply specifically to the thermometer actually 
used, consisting of a film of Aquadag on a thin mica substrate, but the 
same general behaviour may be expected for any carbon thermometer. 
An increase in frequency or a decrease in temperature tends to reduce 
the efficiency and to increase the phase lag. 


§ 9. CONCLUSIONS 
(a) The Mean Free Path 


It has been shown in § 6 that only at temperatures above 0-:5°K is the 
mean free path small compared with the size of the apparatus. Because 
of the trend of mean free path with temperature (table 3) it is thought 
that the continuum theory is probably valid at 0-58°xk and is not far from 
the truth at 0-54°K. On this basis, the mean free paths are taken to be 
0-08 cm at 0-54°K and 0-04 cm at 0-58°K. It has been shown in § 7 that 
at temperatures of 0-4°K and below the mean free path is long compared 
with the tube diameter, and that about one-third of the phonons are 
scattered diffusely at the walls, the remaining two-thirds suffering 
specular reflection. 

Landau and Khalatnikov (1949) have calculated the expected mean 
free path / of phonons, and for temperatures below 0-8°K they find 


J= {1-1 x 1087"9/2 exp (— 4/RT)+5-6 x 10779}! cm, 5 ED) 
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where A/k=9-6 deg. according to Landau (1941, 1947). This function 
together with the experimental points at 0-54°K and 0-58°K is illustrated 
in fig. 11. The experimental points are below the curve by a factor of 
about 6, whereas it is thought that the experimental error is certainly 
not more than a factor of 2. The experimental points of Whitworth 
(1955), from measurements of steady heat flow. and of Kramers e¢ al. 
(1954) from measurements on heat pulses are also shown. Neither value 
is claimed to be accurate, the uncertainty in Whitworth’s results being 
explicitly quoted as a factor of about 4. It should also be remarked 
that the numerical constants in eqn. (19) must be treated with some 
reserve since they depend rather sensitively on experimental data which 
are not very accurately known. 

Fairbank and Wilks (1954) have also measured steady heat flow in a 
tube and they find that their measurements are consistent with a phonon 
mean free path considerably greater than 0-29 mm, the diameter of their 
German silver tube, at temperatures below 0-6°K. They find an effective 
mean free path for thermal conduction of about 0-3 mm at 0-6°K and 
about 0-4 mm at 0-26°K and they remark that this behaviour suggests an 
increasing amount of specular reflection as the temperature is reduced. 
Using eqn. (15) we find from their data €=1 at 0-6°K and €=0-87 at 
0-26°K. Since a quartz tube has been used in the present work and a 
metal tube was used by Fairbank and Wilks, quantitative agreement 
about the degree of diffuse reflection is not to be expected, but the trend 
with temperature is the same in both cases. 


(b) Relation to Other Experiments 


Previous work on second sound below 1°K has all been carried out by 
the pulse method. It might be possible, using the results of the present 
experiments, to predict what pulse shape ought to be obtained, but this 
would involve a knowledge of the frequency characteristics of the 
thermometer and amplifiers used in the pulse work. It is, however, 
possible to compare the values of c’ given in fig. 7 with the velocities 
derived from pulse experiments by measuring the time of arrival of the 
beginning of the pulse. The values at 2 ke/s, for example, are in good 
agreement with the pulse values of de Klerk e¢ al. (1953, 1954) at all six 
temperatures. It seems possible that although higher frequencies were 
present in their pulses, the adverse effect of falling thermometer sensitivity 
with increasing frequency (see §8) prevented frequencies higher than 
about 2 ke/s from being observed. If in fact c’ increases with frequency, 
the front edge of the received pulse ought to arrive at a time corresponding 
to the value of c’ appropriate to the highest frequency detectable. If 
this view be correct then at temperatures lower than 0-5°K an ideal 
thermometer would show a pulse with a front edge occurring at a time 
corresponding to c, the velocity of ordinary sound. The velocities given 
by Kramers et al. (1954) are larger than those of de Klerk et al., a 
fact which might be ascribed to a better frequency response of their 
thermometer. 
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It seems important that hitherto no interpretation of pulse shapes has 
taken account of reflection of the waves at the ends of the tube. The 
present observations are consistent with a reflection coefficient of unity, 
and the amplitudes observed in the shortest tube can hardly be explained 
on any other basis. If, therefore, the leading edge of a pulse arrives at 
the thermometer after a time t, the reflection of this edge will be expected 
to return there after a further time 2¢. Caution must therefore be exercised 
in interpreting the shape of that part of a received pulse which arrives 
after a time 3¢ from the original transmission. 
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XXIX. Collective Effects in Nuclei of Mass 18 and 19 


By F. C. BarKer 
Atomic Energy Research Establishment, Harwellt 


[Received November 10, 1955] 


In an intermediate coupling shell model calculation for the nuclei of 
‘mass 18 and 19, Elliott and Flowers (1955) determined the positions of 
the low lying levels in the absence of any collective effects. They found, 
however that certain E2 transition rates in 170 and 19F could not be 
explained without collective motion and that the introduction of weak 
coupling of the Bohr and Mottelson (1953) type between the surface 
and the nucleons outside the 1°O core gave simultaneous agreement 
in 70 and }°F. In this note we consider the effects of such weak coupling, 
and show that although it enhances the E2 transition rates in 19F by a 
factor of order 5, the changes in level positions and magnetic moment: 
are negligible. 
The condition for validity of the weak coupling approach is 


: 5 k 
AG la ul (es) 4/ (hal) su 


(cf. BM, eqn. (I1.14)). With BM’s hydrodynamical estimates of C and iw 
(about 30 and 7 Mey respectively), EF fitted the lifetime of the 0-872 Mev 
state of 170, which decays by E2 radiation, using 7\/70-2. This implies 
such weak coupling that it is valid to consider only states containing 0 or 
1 phonon, and to use perturbation theory to calculate changes in wave 
functions and energy values. 

To fit the quadrupole moment of the ground state of 170, a much 
smaller value of x\/j is required. This could be attributed to the matrix 
elements of k(r) between different single particle states being different, 
contrary to the assumption of BM. This may be due to the single particle 
states involved being relatively loosely bound (a different reason has 
been suggested by Blin-Stoyle (private communication)). In order to 
fit the experimental values for the 17O ground state quadrupole moment 
and 0-872 Mev state lifetime 


Q=—0-005 x 10-** cm?, =. r= 2:5 X 10719 sec, 


with the nuclear radius R,—1-44A13 10-3 cm, we require (cf. BM 
eqn. (II.10)) 
(id |k|1d)=0-033C, (1d | k | 28)=0-200 
eS Le oe Ue 
+ Communicated by B. H. Flowers, 
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(the matrix elements are assumed to have the same sign). This coupling 

alters the positions of the levels of 170, and EF’s values for the strength 

of the spin-orbit force and the d-s level separation are replaced by 
é=2 mev+0-0010C0, d-s=1-125 mev—0-0052C. 

For the mass 19 nuclei, the unperturbed wave functions are taken 
from those determined by EF, with only the maximum orbital symmetry 
L-S states retained for 7’=4, and only the (d;/.3~*s”) j—j states retained 
for 7 =%. The relative changes in the energies EH, of the lowest states 


with given 7' and J are found to be 


A(E y/o, 5/2 1/2, 1/2)=0-0016C, A(Ex;2, 3/2 —E a2, 1/2) =0-00180, 
A(E3;9. 5/2 1/2, 1/2)=0-0065C, A(Es;, 3/2 —Lyy2, 1/2) =0-00430, 
A(E se, v2—Ly2, 1/2) =9-0070C, 


and the change in the 1°F ground state magnetic moment (in units of 
nuclear magnetons) is 
C 


Ap=—0-0011 —. 


For reasonable values of C and hw, these changes are negligible. 
Recently the magnetic moment of the 0-197 Mev 8+ state of 1°F has 
been measured. The results are 


p=4:0+0-9 (preliminary value, Lehmann et al. 1955), 
o=3:4+0-9 (Treacy, private communication). 


EF’s wave function for this state gives u.=3-3. The change due to 
surface coupling, 
C 
Au=—0 0016 5— ; 
is again negligible. 

K2 transitions in !°F are enhanced by the surface coupling by a factor 
of order 5, but those in competition with M1 transitions between the 
same states are still unimportant. Even in the 3++4+ transition for 
which EF found an M1 transition rate about 1% of the single particle 
value, the E2 transition rate is only about one-fifth of this. * 
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XXX. The Lowering of Fracture-Stress due to Surface Adsorption 


By N. J. Pretcu 
Metallurgy Laboratory, University of Leedst 


[Received September 15, 1955] 


ABSTRACT 


Gibbs’ adsorption equation and the Griffith and dislocation theories 
of fracture are used to calculate this lowering. The possibility that such 
adsorption causes the embrittlement of iron by hydrogen is discussed and 
is shown to be in good agreement with measured fracture stresses. 


§ 1. InrRopUCTION 

Tue interest in this topic dates from the suggestion (Orowan 1944) that 
adsorption of atmospheric constituents on surface Griffith cracks in 
glass lowers the fracture stress and produces the delayed-fracture effect. 

The stress o, required to enlarge a surface crack of depth ¢ in a solid 
of surface energy y and Young’s modulus ZF is 

o7=(2Hy/ne)'?. ee ieee pete ee Ly) 

If the surface of the crack is covered by an adsorbed phase, y is lowered 
to y, and o, is correspondingly lowered, provided the extension of the 
crack is sufficiently slow to allow the formation ae increments of surface 
that contain the adsorbed phase. 

Thus, a knowledge of the surface energy lowering is required for 
calculation of the new fracture stress. 


§ 2. THe Surrace Enercy LOWERING 
This is, of course, given directly by Gibbs’ adsorption equation 
dy=—TI dp, 
where dy is the change in surface energy produced by a change dy in 
the chemical potential of the adsorbed phase and I’ is the number of 


molecules adsorbed per unit area. 
If the adsorption takes place at temperature 7’ from a gas at pressure /, 


, dy=—IkT din p. Rane ah Deere the (2) 


From the Langmuir isotherm, 
Ap 
1A py ; 
where I’, is the number of molecules adsorbed per unit area at saturation 
and A is a constant. 


+ Communicated by the Author. 
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With a diatomic molecule that dissociates on adsorption, I, is halved 
and the isotherm becomes 
(Ap)! 
1+ (Ap) 
Substitution for I’ in (2) and integration gives 
yi=y—UkT In (1+-Ap), 
or, for a diatomic molecule that dissociates, 
yy=y—2T kT In (1+(Ap)”?). 5 Ne) ioe een oe 
The lowered fracture stress can then be obtained from (1). 
The stronger the adsorption, the greater is A and, as would be expected, 
the greater is the fracture-stress lowering. 
If the Langmuir isotherm is not sufficiently accurate, it may be 


necessary to use one of the other isotherms or a numerical integration 
of the experimental one in this calculation. 


§ 3. THe FRAcTURE OF METALS 


In this case, we are fairly certainly concerned with the production of 
cracks ahead of arrays of dislocations held up at some obstacle, rather 
than with the enlargement of pre-existing Griffith cracks. 

Stroh’s (1954) criterion for the production of a crack can be put in the 
form 

T=12y/nb, 


where the shear stress t compresses an array of n edge dislocations of 
Burgers vector ). 
From Eshelby, Frank and Nabarro (1951), 


n=n7(1—v)rL/|pb, 


where L is the length of the array and w is the rigidity modulus. 
Thus, the tensile stress required at 45° to the plane and length of the 
array to produce a crack is 


‘ l2wy i 
St a ee : Maven ae, Soom 
“ (= fe . 


With real metals, a considerable stress oy may be necessary merely to 
move a dislocation in the crystal, so og must be added to o in (4) (Hall 
1951, Petch 1953, Cracknell and Petch 1955 a). This equation then 
gives the criterion for the cleavage fracture of polycrystalline specimens, 
and it seems (Petch 1955) that it will also apply to ductile fracture, at 
least for «-iron. Thus, the fracture stress is 


#: Byy i). 
o,=0,-7-4( ———_ : 
four (<5) . vets Bee Oe 


It will be appreciated that this calculation assumes isotropy. 
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If a very mobile solute is present, so that the cracks can form with 
adsorbed solute on them, ' 


3 1/2 
oy=o9+4 (= {y—DkT In (1+-4p)}) 


or, for a diatomic molecule that dissociates on adsorption, 


3 1/2 
te aes ( (y—22kT In (1+(4p)"™)}) ? - (6) 


where p is the equilibrium pressure of the solute. 

For polycrystalline «-iron and zinc, the experimental evidence indicates 
Lol, the grain diameter (Petch 1953, Greenwood and Quarrell 1954). 
In fractures at small plastic deformations, L is probably 1/2, and, in the 
case of «-iron, it appears that the effective value of Z is not altered by 
plastic deformation (Petch 1955). 

From (6), with Locl, the lowering of (c;—o9) produced by the 
adsorption depends upon the grain size, being proportional to /-1/2. 


§ 4. HyDROGEN-EMBRITTLEMENT OF [RON 


A small quantity of hydrogen in iron markedly lowers the fracture 
stress. It is quite widely believed (Zapffe and Sims 1941) that the 
accumulation of hydrogen in small cavities present as defects is the 
cause ; it is thought that the stress around the cavities initiates premature 
fracture. The room-temperature equilibrium pressure of hydrogen is 
certainly very high, even when only small amounts are in solution, and 
it may well be that these internal stresses do explain the suppression of 
the yield point that hydrogen can produce (Cracknell and Petch 1955 b), 
but their role in fracture is less certain. The plastic deformation prior 
to fracture will tend to homogenize the shear-stress distribution in the 
metal. 

Since an iron fracture surface will be a very strong chemisorber of 
hydrogen and since hydrogen can diffuse extremely rapidly through iron, 
Petch and Stables (1952) suggested that adsorption might be the cause 
of the reduced strength. In keeping with this, hydrogen has no effect 
on the strength at high strain rates. 

The movement of hydrogen involved in the production of a crack at 
an array should be roughly comparable to the processes of quench- and 
strain-ageing, so, from the known diffusion rates, it should require 
~10-5 sec at room temperature. This makes it a reasonable possibility 
in a normal tensile test. 

Equation (6) can now be used for a quantitative evaluation of the 
effect that hydrogen can produce. 

The pressure of hydrogen (mm of Hg) in equilibrium with «-iron 
containing scm? of hydrogen/100 g is given (Sykes, Burton and Gegg 
1947) by 

pi/2—0-36s exp (3650/7). 


Z2 
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One way of obtaining A in (6) is by noting that it equals p; fe oa yess: 
P1/2 is the hydrogen pressure for [/[',=}. This gives A~5x 104 mm— 
for hydrogen on nickel at room temperature, and the same value was 
used for hydrogen on iron, since the two adsorptions are very similar, 
with practically the same heat of adsorption (Beeck 1940, 1 950). 

Figure 1 shows the lowering of the surface energy of iron at room 
temperature calculated from (6) using this data. The effect is large. 


Fig. 1 


ERGS/SQ. CM. 


SURFACE ENERGY CHANGE 


c.c.H,/ 100 cm 


Lowering of the surface energy of iron by adsorption of hydrogen. 


To obtain the corresponding values of (c;—o)) when hydrogen is 
present, the initial surface energy is required. A value of 1180 erg/em? 
was obtained by the application of (5) to previous measurements (Petch 
1955) of the fracture stress/grain size relationship for a hydrogen-free mild 
steel, assuming L=1/2, Taking a medium grain size of [-¥/2=5 mm-1/2 
as an example, a fall of (c;—o 9) from 20-5 tons/in.2 in the absence of 
hydrogen to 9-5 tons/in.? at 10 em? H,/100 g is then indicated. 

The influence of temperature on the effect of hydrogen adsorption can 
be estimated by using the Clausius—Clapyron equation to obtain the 
adsorption isotherms for other temperatures from the room temperature. 
one. The heat of adsorption is somewhat uncertain. It falls from an 
initial 33 kilocal/g mole as the fractional coverage @ of the surface 
increases ; 15 kilocal is reached at 0-1 mm and room temperature, which 
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probably corresponds to 6~0-75, and presumably the value at 6~1 is 
only a few kilocal (Beeck 1950). A calculation for —115°o using this 
information suggests that the change from room temperature does not 
have much effect on the surface energy lowering. 


$5, COMPARISON WITH EXPERIMENTS 

Figure 2(a) shows the dependence of the fracture stress of a poly- 
crystalline mild steel containing ~10 cm? of hydrogen/100 g upon the 
grain diameter / at 18°C. Figure 2(b) gives the previously-determined 
figures for the same steel without hydrogen (Petch 1955) 


O; TONS /SQ.IN. 
KG./SQ. MM, 


V2 —Y2 
tL MM. 


Fracture stress/grain size relationship for a mild steel at 18°c, 
(a) with ~ 10 cm? hydrogen/100 g, (6) without hydrogen. 


In both cases, the relationship is of the form 


where o, and k are constants. It is apparent that hydrogen causes a 
reduction in the slope k and in a». 

The reduction in k& is consistent with a lowering of the surface enerey, 
as in (6). AtJ-1/2—5 mm~1?, the observed value of (o,—a») is 11 kee 
in the presence of the hydrogen, compared with the figure of 9-5 tons/in. 
obtained using the calculated surface energy lowering. aie 

The reduction in co, by the hydrogen is readily explained as a direc 
result of the lower strain at fracture, which, in turn, 1s caused. by the 
reduction in &. The op term simply represents the stress required to 
move a dislocation and this movement will be easier if the strain-hardening 
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is reduced. Thus, oo will fall as the average reduction of area falls from 
67°%, to 21% in the presence of the hydrogen. 

Figure 3 shows that lowering the temperature to —115°c has little 
influence on the hydrogen effect. The average reduction of area is 
reduced from 61°% to 23°% by the hydrogen. 


Fig. 3 


KG./SQ.MM. 


O, TONS/SQ.IN. 


Fracture stress/grain size relationship for a mild steel at —115°c, 
(a) with ~ 10 em’ hydrogen/100 g, (b) without hydrogen. 


All these measurements are in good agreement with the explanation of 
hydrogen embrittlement in terms of adsorption. 


§ 6. EXPERIMENTAL DETAILS 

The general procedure followed that used in the measurements on the 
same mild steel when free from hydrogen (Petch 1955). 

In the present case, the tensile specimens were machined to 0-15 in. 
diameter on the gauge length and annealed at 650°c. The hydrogen 
was then introduced electrolytically using a 4% H,SO, electrolyte 
containing a few drops of yellow phosphorus in CS, as a poison. The 
current density was 2 amp/in.?, and it was applied for 17 hours. Since 
electrolytic charging may produce surface blisters, the specimens were 
turned down to 0-125 in. immediately after charging; 0-001 in. cuts 
were used and the specimens were oil-cooled during the machining. 
They were subsequently polished to a 3/0 emery finish and fractured. 
Some specimens were machined directly to 0-125in. gauge-length 
diameter, annealed, hydrogen-charged, polished and tested ; the fracture 
stresses agreed with the previous ones. 


! 


Lowering of Fracture-Stress due to Surface Adsorption 337 


Within each series of measurements, there was some variation in the 
strain at fracture with grain size. Since this introduces a variation 
i go, a small correction was necessary to bring the results to an average 
strain. To obtain this correction, cg was plotted against /-1/2 without 
any correction, and the difference in 0) between the hydrogen-containing 


and the hydrogen-free specimens was then used to find the variation 
of o, with strain. 
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XXXI. The S-matrix for Neutral PS-PV Meson—Nucleon Interaction 


By J. 8S. R. CHISHOLM 
University College, Cardiff + 


[Received October 28, 1955] 


SUMMARY 
Nucleon and neutral meson fields with an exponential interaction, 
derived from the PS—PV interaction by the Dyson transformation, are 
studied. A new system of graphs called ‘skeletons’ is described ; the 
S-matrix for any process is the sum of terms corresponding to skeletons 
with the appropriate external lines. 


$1. IyrropucTION 

i. J. GLAUBER (1951) has investigated neutral meson theories in which 
an exponential of the meson field occurs in the interaction Hamiltonian. 
This paper generalizes Glauber’s results by deriving the S-matrix for an 
exponential interaction with fermion and neutral meson fields both 
quantized. The meson field considered is pseudoscalar, and the inter- 
action Lagrangian is derived from the neutral PS—PV Lagrangian by the 
Dyson transformation (1948); the transformation is given in § 2, and 
involves a mass renormalization. In §§3 and 4 a graphical theory of 
the interaction is developed, and in § 5 it is shown that the renormalized 
S-matrix is a sum of terms corresponding to a new system of graphs called 
‘skeletons’. The fermion and meson propagators in the skeletons 
contain sinh 4,(x) and cosh 4,(2), and correspond to the propagation 
of an arbitrarily large number of mesons. Section 6 consists of a discussion 
of a further renormalization required to eliminate fermion lines which 
begin and end at the same vertex. Section 7 consists of a brief discussion 
of possible interactions. 


§ 2. THE TRANSFORMATION OF THE LAGRANGIAN 
The Lagrangian density for neutral meson and nucleon fields with 
PS-PV interaction is 
L=L,4+L1,+2L,+1,; 
where L1= (a) y,0,3)(@), 
Ly=Mip(x)p(x) 
L3=—31(0,9)? +m $"], 
Ly=igh(x)y5ya(%)(O,P)s 


with p(a)=p*(x)yy. 
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We perform a unitary transformation (Dyson 1948) on the system — 
by putting 
aan. [tGy sb(a)|f(a) 4 
and P(x (x) exp [—iG'y54(z) 
the y,; operating on the spinors. 
This transformation leaves 1, unchanged, but affects on L, and L,: 


L,=6*(x) exp [—iGy,4(x) lya(y,0,,) exp [iGy h(a 


es 


(1) 


Now 
d, exp [tGy,¢]= aaa ae (Gear = 


ee exp iG byt (0,¢6)—4G°[4, 0,0 | exp [(Gdys] ; 
the quantity [4, 0,4] is infinite for .=4, as it contains a $-function with 
zero argument. We have 


L,=C* (x) yay,diGy (0,0) — bee Bee C(x 
+Tlay,atla); tietaet ea0 (2) 

also 

Ly=190* (x) yay sy 6 (%)(0,8); 
which is cancelled by the first term in (2) if G=g. The second term in 
(2) is an infinite mass renormalization term of the type discussed by 
R. J. N. Phillips (1954), and is dropped, as is usual for renormalization 
terms. Then L, and L, reduce to 


Ly =L(a)y,0,£(20). Dee: i ee (3) 
Now 
L,=ME*(x) exp [—iGy sly, exp [tGy56]¢(x) 
=L,'+ME(x) {exp [2igy 36(x)]—1 }<(a) 
where LM C(x) C(x). 
Thus L=L,4+L,/+L,'+L, 
where DL, =MU(x)[exp (21gy5¢)—1] d(x) ; Oe Me), 


L,{ is the transformed interaction Lagrangian. 


§ 3. GRAPHICAL REPRESENTATION 


The new interaction term (4 a) is expressible as a power series in ¢(2) 


Ly =MU(u) E (2igysh(e)]"C@)in | 
n=1 
=MU(w)dle) ¥ Pigdee P42! 


+ MC(a)y 5h(w oa B [2igd(x)]* 1/(21—1) !. tent D) 


We can use the familiar caer mee of perturbation field ee 
to treat the series of perturbation terms (46). A term containing [ d(x) 
corresponds to a vertex with an incoming and an outgoing ae 
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line, representing {(x) and {(x), plus n meson lines. These vertices fall 
into two classes, ‘ vertices of odd parity’ with n=2/—1, in which a 
factor y, occurs at a vertex, and ‘ vertices of even parity’ with n= 21, 
in which no such factor occurs. Pairs of operators 2ig¢(x,), 20gP(x;), 
and C(x), t(a;) from vertices at x; and w,; can be associated together in 
P-brackets, giving rise to multiples —4g?4,(a;, #;) and S,(x;, 2;) of the 
usual Feynman propagators, and we can draw graphs in which vertices 
are connected by internal meson and fermion lines, these lines being 
associated with the appropriate propagators. In addition there occur 
vertex factors M(y;)"/n!. Clearly two vertices with a large number 
of meson lines emanating from each may well be joined by many meson 
lines, and we must also remember that meson lines may begin and end 
at the same vertex. As in the simpler theories with scalar interactions, 
the S-matrix element for any process is the sum of terms corresponding 
to all graphs with the appropriate external lines; but the system of 
graphs is very complex as any number of meson lines may connect each 
pair of vertices. In order to classify the system of graphs we introduce 
the concept of the ‘skeleton’ of a graph. Consider any graph; in 
general there will be more than one line (meson or nucleon) connecting 
any two vertices. The ‘ skeleton’ of the graph considered is formed by 
rubbing out internal meson lines until only one line (meson or nucleon) 
joins each pair of vertices that were joined originally. The set of all 
graphs with the same skeleton will often be considered together in the 
following work ; two or more graphs with the same skeleton will be called 
‘similar graphs ’. 


§ 4. ENUMERATION OF GRAPHS 

Sorresponding to each graph there will generally be many similar 
terms in the S-matrix, since there are many ways of selecting the pairs 
$(x,), $(x;) from the nth and pth terms [¢(x,)]" and [(x,)}? of Ly’(x,) 
and L,'(x;). We shall now evaluate the weight of any graph, that is, 
the number of terms in the S-matrix which correspond to it. 

Consider a vertex in any graph with » meson lines emanating from it. 
Let these lines be grouped as follows : 

(a) a group of 2g lines which join up in pairs to form ‘loop lines’ at 
the vertex ; 


(6) a group of n, lines which are external meson lines ; 
(c) (m—I1) groups of n, lines (r=2,...,m) which proceed in these 
groups to (m—1) other vertices. 
m 
Then 2g+ 2 n,=n. 
r=1 


Regarding each line within a group as indistinguishable, the number 
of ways of choosing the grouping of the n lines is 


n! 


(2q) 1 TIn,! . . : . . . . . (5) 


Neutral PS-PV Meson—Nucleon Interaction 341 


The number of ways of joining up n, internal lines which leave one 


vertex and mn, lines which leave another, or of arranging n, external 
lines, is 


N, 4 ee he Satse oI 4.4 (G) 

the number of ways of joining up 2g lines in pairs is 
(29) ! 
Dat ees a ay eer ne (7) 


The product of a factor (5) for each vertex, a factor (6) for each 
group of nm, meson lines (r=1,...,m), a factor (7) for each vertex, 
and a factor (n!)~! for each vertex from (46), gives rise to the 
following product of factors: (i) (n,!)~! for each group of n, internal 
meson lines ; (ii) (2%)~1 for a vertex with q loop lines. 

There are also in the matrix element factors (iii) M(y;)" for a vertex 
with m meson lines ; (iv) [—4g?4,(z,, x;)]"" for each group of », meson 
_lines joining vertices at x; and x,; (v) [—4g?4p(x,, x,)|* for ¢ loop lines 
at a vertex ; (vi) Sp(x;, x;) for a fermion line joining vertices at x, and 2; ; 
(vii) external line factors of the form C(x), C(x) and (x) ; (viii) a factor 
(—1)*, where & is the sum of the number of vertices in a graph and the 
number of polygonal fermion ares (open or closed) ; this factor is due to 
anticommutation of the fermion operators, and is the same for all graphs 
with the same skeleton. . 


§ 5. Murtti-Mreson PRoPpAaGATORS 
We are now in a position to sum over S-matrix terms corresponding 
to similar graphs ; this is equivalent to summing over all possible n, for 
each group of internal lines. For the moment we ignore the factors (ii) 
and (v), which merely give rise to a renormalization. At first sight the 
factors (i) and (iv) above give rise to meson propagators something like 
exp [—4g?4,(a;, x;)]. However the factors (iii) depend on 


N=N1 +. - tN +2¢ 5 
that is, they depend on the number of meson lines in each group at the 
vertex. It is therefore necessary to distinguish between terms in the 
meson propagators with x, even and those with n, odd. We define odd 
and even meson propagators d,(v;, z;) and d,(a;, ;) by summing over odd 
and even n, respectively : 


A(X ;, %;)= 2’ [—49?4 p(x, a) 4(21-— 1) ! 
t=1 


=sinh [—4g?4 p(x,, x;)], ee ee So) 
A(X; Xj)= eee [4922 (x; @,)"} (21) | 
=e cosh: (= 497A pe ;)| ee ee ee (OU) 


likewise we define odd and even fermion propagators 
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8,(x;, 2;)=sinh [—4g?4p,(x;, 2;)] . Sp(x;, &;) seen (8 a} 
and 8,(%;, %;)=cosh [—4g?Ap(x;, X;)].Sp(w;,%;). - . . (9 b) 
Skeletons can be drawn in which any number of fermion lines are joined 
end to end. - To take account of this iteration of fermion lines, we should 
replace s, and s, by ‘iterated propagators’ S, and S,, defined by a 
Bethe-Salpeter equation. In momentum representation, S, and S, are 
simple algebraic functions of s, and s,. However, the Fourier transforms 
of s, (w) and s, (x) are infinite ; this difficulty must be overcome before 
S, and S, can be properly defined. 

Now take all terms formed from the skeleton by writing (a) d, or d, for 
each meson line; (b) s, or s, for each fermion line; (c) M or My, for 
each vertex (according to the parity of the vertex, as discussed below) ; 
(d) ZT, C and ¢ for the appropriate external lines ; then the sum of all such 
terms is the sum over all graphs with the same skeleton, apart from the 
factor (—1)*, as it takes into account all possible values of the , for 
each group of meson lines. It remains for us to specify the factor (c) 


m 

more clearly ; a factor y; occurs at a vertex when 2 vn, at the vertex 
r=1 

is odd. If we are evaluating a particular term with various odd and 


even propagators corresponding to lines in the skeleton, it is best to 
label the lines with minus (odd) and plus (even) signs, giving each a 
parity. - The parity of a vertex in the skeleton is then the product of the 
parities of the lines arriving at the vertex ; then it is easy to insert the 
correct factors My, and M for the odd and even parity vertices 
respectively. 

The factors (ii) and (v) have so far been left out of account. They can 
be treated separately, since the numbers q occur only in these factors, 
(y;)*4 in factor (iii) being equal to unity always. These factors, summed 
over q at a vertex, give a term 


ive) 


Eg! —892A p(x;, &,) exp [—892Ap(a,, 2]... (10) 


q=0 
Since 4,(v, x) possesses a 5-function infinity, the vertex factor (10) is a 
zero of a high order. It is however just a renormalization factor which 
appears at every vertex. Without any more justification than is provided 
by similar situations in simpler theories, we omit these vertex factors. 


§ 6. FermMIon Loop Lings 
A further type of divergence is encountered when we consider vertices 
in skeletons which have several meson lines (either external or connecting 
with other vertices) and a fermion ‘ loop line ’ which begins and ends at 
the vertex. The cloud of mesons which usually accompanies the fermion 
line has been eliminated already by the meson loop line renormalization : 
so the fermion loop line contributes a factor } 
dp 
typ+M 
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to the matrix element ; this is a divergent integral which is also en- 
countered as a vacuum polarization term. It is independent of the 
number and the nature of the meson lines at the vertex, and it can 
be subtracted out completely by introducing a term 


— an te | Set | _4p ] 
n! iyp+M 


into the Lagrangian, where n is the number of meson lines at the vertex 
in the graph considered. Summing over n, we see that all skeletons 


containing fermion loop lines are entirely eliminated by adding to the 
Lagrangian a term 


== 'Tr [exp (2igy 50) | tl ; 


This further renormalization enables us to omit skeletons containing 
fermion loop lines. 


§ 7. EXAMPLES 

First we shall demonstrate by a simple example the type of graph and 
skeleton that can occur. The example is what we might term a third- 
order nucleon-nucleon interaction. The graph (a) represents the type 
of graph we are considering, with a cloud of mesons passing between each. 
pair of vertices. The aoa of all such graphs is shown in diagram (5). 
the RNA being labelled 1, 2, 3, for convenience. Each line can have 
parity + or —, so that ee are eight possible labellings of the line 
parities ; we have labelled diagram (>) with one of the eight labellings. 
Diagram (c) shows the corresponding vertex parities. From the skeleton 


(a) Cb) (ee) 


{ bu uy 
my 4 


@ | eae eae 
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we can write down the matrix element ; it is a multiple of 
[[ [dary dra desl Elen) Ms,(ey, 24) My oL(e1)] 


x [C(ag) My 56(a3)] d-(%2, 3) Ao(23, 73). fy 20" cere nee 


The questions of convergence and evaluation of integrals such as (11) are 
complex, and no attempt is made here to discuss them. 

As in electrodynamics and other weak-coupling theories the simplest 
processes (represented by the simplest skeletons) do not involve evaluation 
of integrals in momentum space of the type which occur in (11). The 
simplest process we can consider is meson—nucleon scattering in first 
order, with skeleton (d); we do not even have to deal with the new 
propagators to work out this term. The simplest nucleon-nucleon 
scattering process has skeleton (e) ; at low energies we may expect only 
the even propagator to contribute to this term, because of the y, factors 
occurring with the odd propagator. The simplest meson production 
skeleton (f) is very similar to (e); skeletons of this type also give rise to 
multiple meson production. Glauber (1951) has already discussed 
physical effects corresponding to the simplest graphs. 


§ 8. CONCLUSION 
This method of dealing with gradient interaction is not readily extended 
to charged and symmetric meson theories ; the Dyson transformation 
does not apply to these theories, due to the non-commutative properties 
of the r-matrices. It may be possible however to define a new symmetric 
theory by generalizing the interaction Lagrangian (4 a). 


Since this paper was submitted for publication, some work by 
R. Arnowitt and 8. Deser (1955) has come to the notice of the author. 
This work covers much the same ground as the present paper, and 
the two pieces of work have been done simultaneously. 

The author is very grateful for much help and constructive criticism 
given by Dr. K. Burton, Dr. A. H. de Borde and Dr. R. G. Moorhouse, 
and for the constant encouragement of Professor J. L. B. Cooper. 
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SUMMARY 


The flow of a viscous fluid past a sphere at high Reynolds number as 
given by Oseen’s approximation, in which the inertia terms of the 
equations of motion are linearized, is discussed. It is found that a 
boundary layer is formed on the forward side and a wake extending to 
infinity on the rearward side. Inside the wake, which is bounded by a 
circular cylinder with its generators parallel to the undisturbed direction 
of the streamlines, the fluid at infinity is moving towards the sphere, 
while outside it is moving away with the undisturbed velocity, and there 
is a shear layer at the cylinder. 


§1. InTRODUCTION 

RECENTLY in two papers (Pearcey and McHugh 1955, Tomotika and Aoi 
1950) the steady flow past a sphere at low Reynolds number R according 
to Oseen’s approximation has been discussed. In the first of these 
Tomotika and Aoi (1950) consider the flow at R=2aU/v=1, where U is 
the velocity of the undisturbed stream, a is the radius of the sphere and 
v is the kinematic viscosity, and announce the presence of a small vortex 
behind the sphere. Their calculations are extended in the second paper, 
by Pearcey and McHugh (1955), who however do not find such a vortex 
even when R=10. There is thus some doubt whether Oseen’s equations 
of motion will reveal the presence of a vortex at any Reynolds number. 
Since a solution of these approximate equations is useful as a guide to the 
solution of the full Navier Stokes equations it is of interest to examine it 
at high Reynolds number to see if the vortex is present. 

A full discussion of Oseen’s equations at high Reynolds number was 
first given by Oseen himself (1927) who obtained the conditions satisfied 
by the inviscid part of the flow in general and obtained the explicit 
solution for the steady flow past a circular cylinder for which there is such 
a vortex and a thin plate. His discussion is in some ways a little involved 
and so here in § 2 a shorter derivation of his principal results is given. In 
§3 these are applied to the sphere in a uniform flow. 

Let C be a circular cylinder circumscribing the sphere and extending 
downstream of it, and with its generators parallel to the x-axis which is 


+ Communicated by the Author. 
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the direction of the stream. It is then found that inside C the fluid is 
moving with vorticity, and with a speed of never more than about a quarter 
of the mainstream velocity, towards the sphere and C. It reaches C first 
and there the velocity component parallel to the x-direction is discon- 
tinuous so that outside C the fluid moves away from the sphere with a 
velocity which is ultimately constant. Beyond this fluid moving in an 
infinitely long vortex is fluid which has come directly from an infinite 
distance upstream and whose motion is irrotational except in the im- 
mediate neighbourhood of the forward side of the sphere where there is a 
boundary layer. The thickness of the boundary layer is O(v) as against 
O(v/?) for a boundary layer derived from the Navier Stokes equations. 
Of course the discontinuity on C is the limit of a shear layer whose thick- 
ness is O(v/2) and it may be regarded as an extension of the boundary 
layer on the forward side of the sphere. The presence of the boundary 
layer on the forward side of the sphere confirms the trend indicated by 
_ Pearcey and McHugh (1955) but they did not find the vortex and from 
their diagrams there were only the barest of hints that it would ultimately 
appear. 

It is interesting to consider whether the approximate solution at large 
values of R throws any light on the exact*although unstable) solution of 
the full Navier Stokes equations. For this purpose we must bear in mind 
first that the basic assumption of the approximate theory is that any 
transportable quantity such as vorticity is carried not with the fluid but 
parallel to the w-axis with the undisturbed velocity of the fluid. Thus it 
is seriously in error in the boundary layer on the forward side of the sphere 
leading to the wrong order of magnitude for its thickness, and in the 
wake where the direction of motion of the fluid is almost reversed. Second 
the approximate equations being linear there is no breakdown in the 
<lerived boundary layer equations when the skin friction vanishes and in 
fact the only restriction on the boundary layer is that it is confined to the 
forward side of the sphere. For a fluid satisfying the Navier Stokes 
equations a paper of Goldstein (1948) indicates that when the skin friction 
vanishes the conventional boundary layer probably comes to an end and 
there is ample evidence from observations to confirm it. The main 
stream then leaves the body. With these two differences in mind it is 
inferred from the approximate solution that near the front stagnation 
point there is a boundary layer which separates from the sphere when 
the skin friction vanishes. Outside the stream tube which leaves the 
sphere here is the main body of the fluid moving irrotationally ; inside 
the flow is vortex like and may satisfy the no slip condition on the sphere 
without any necessity for a boundary layer. 

The simplest solution of this type is that in which the fluid inside the 
breakaway stream tube is at rest (Kirchhoff flow) and for a circular 
cylinder the corresponding flow has already been suggested as being 
appropriate (Imai 1953, Kawaguti 1953, Squire 1934). However such a 
conclusion has a grave objection for the two dimensional flow. The 
argument in favour is briefly as follows. Let us suppose that breakaway 
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occurs at an angular distance « from the front stagnation point. Then 
Schmieden (1932) has shown that if «<55° there is a contradiction because 
the free streamline crosses the cylinder. If «>55° the free streamline 
has an infinite curvature at breakaway which means an infinite pressure 
gradient so that the skin friction vanishes at an angular distance B<a 
from the front stagnation point, which is also a contradiction. Further 
as a> 55°-++, B+ 55°+ and so it is concluded that the correct value of « 
is 55°. The objection to this is that if «55° the pressure gradient up- 
stream of separation is always favourable so that at 55° the skin friction 
is certainly not zero. Although it has not been proved that breakaway 
cannot occur before the skin friction vanishes (so that it would be dis- 
continuous there) no evidence in support is forthcoming and it is unlikely. 
In that case the limiting argument given above is misleading and ir- 
relevant once it is established that the limit flow does not satisfy the 
condition required. 

Another solution of this type which is unacceptable is that given by 
Lighthill (1953) for the flow over a step on a flat plate. He assumed that 
separation occurred at a distance x from the leading edge of the plate and 
calculated the resulting Kirchhoff flow. The value of a was then found by 
an approximate boundary layer theory. The objection to this is that 
whatever value of «x is chosen the pressure gradient is singular at x so that 
separation must occur at a distance y<w from the leading edge. Thus 
if the flow is of the Kirchhoff type separation must occur at or very near 
to the leading edge of the plate. 

A second possible type of flow at infinite Reynolds number may be 
inferred from some recent work of Southwell and Allen (1955) who ob- 
tained a numerical solution of the steady flow past a circular cylinder at 
Reynolds numbers of 1, 10, 100, and 1000. A study of the streamlines 
they display suggests that the flow at infinite Reynolds number may 
consist of a potential flow together with a pair of vortices. This inter- 
pretation of their work leads to a difficulty also, for such a solution would 
involve breakaway at a stagnation point. However in all known instances 
of a main stream accelerating from one stagnation point and then slowing 
down to another the skin friction vanishes well before the rear stagnation 
point. Further numerical work seems called for to elucidate the behaviour 
of the vortices and to ascertain whether perhaps more are formed at 
higher Reynolds number. 


§2. Frow at Large ReyNoLpS NUMBER 

The body is taken to be at rest and at an infinite distance upstream the 
fluid moves in the direction of the #-axis (a increasing) with uniform velo- 
city U. In the form of Oseen’s approximation which we shall use here 
the non-linear terms (u. grad) u, where u is the fluid velocity with com- 
ponents (wu, v, w), in the Navier Stokes equations are replaced by U(du/dx) 

reducing them to 
| du 
Ox 


sER. 8, VOL. I, NO. 4.—APRIL 1956 2A 


ee ad eer divd=s0 5. . veut (254 
P 
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where p is the constant density of the fluid, p is the pressure and v the 
kinematic viscosity. It is noted that Oseen (1927) preferred a slightly 
different form for (2.1) in which p is replaced by g=p+zpu®. The 
difference in physical terms is that (2.1) assumes that the perturbation 
of the velocity is small compared with its undisturbed value whereas. 
Oseen’s form assumes only that the product of the perturbation of the 
velocity and the vorticity is negligible. If the motion is irrotational or 
nearly so it is clear that Oseen’s form is more acceptable and would give 
a more accurate estimation of the pressure. However, such irrotational 
flows are always adjacent to boundary layers in which the velocity 
changes rapidly so that Oseen’s form is bound to lead to pressure dis- 
continuities at infinite Reynolds number which are known not to occur 
in practice. In the form which we choose however the pressure is always 
continuous and so it is felt that a more realistic overall picture is obtained. 
In any case it is only the pressure which is involved in the difference and 
we cannot expect too close an agreement with the correct solution which- 
ever point of view we take. 

The boundary conditions to be satisfied are that at an infinite distance 
upstream 

u=(U,/0;.0) 2 ones 

and on the body u=0. 

In the usual way we write 


p=—pU(ddlox) where V7d=0 .. . « (2.3) 
whence 
us-grad ¢4-ut® +. sats s)he ee eee 
with 
du* 
V7u*t¥=2h—— , ~div u*=0 and k=—U/2y.. . . .. (2.5) 
Ox / 


The principal result we wish to establish in this section is that if 
(%, y, 2) is a point on the body and (a, y, z) a neighbouring point outside it 
u* (x4, Y, 2)=U*(X, y, 2) if wy >2Xp 

and 

u*(z,, ¥,2)=0 2%. = Slee 
in the limit as k>co. Thus on the side of the body facing the oncoming 
stream there is a boundary layer at large values of the Reynolds number 
(i.e. &) in which u* rapidly falls from its value on the body to zero just 
outside, while on the rear side u* is independent of 2 and so a wake is 
formed. These results may be proved fairly easily as follows. Let us 
consider a finite value of #. Then at infinity u*=0 and ¢= Uz to satisfy 
the boundary conditions there. We suppose that the equation of the 
body is =/(y, z) and « represents distance from the body along the out- 
ward drawn normal which has direction cosines (1, m, n) relative to the 
orthogonal frame of reference. If, on the body, 


0 
u*—A(y, 2), an [u* exp (—ka)]=kBly, z)exp(—ka) . . (2.7) 
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then at any external ae 


uta | eS ses (sey BBG By one ees mee) 
a1, 2 [ exp (—kR 
—exp (—ka')A(y’, z 2)a, je | ee 
where R?=(x—wx')?+ (y—y’)?+ (z—z’)? the integral being taken over the 


surface of the body. The Pehaniout of the integral is required when kis large 
it pong assumed and proved a posteriori that A and B do not vary rapidly 
with y’ and 2’ except possibly at isolated points. No evidence for such a 
possibility is present in the solutions at small values of & and as we shall 
see below a complete solution may be found in which A and B do not vary 
rapidly. 
It is sufficient to consider the contribution u,* to u* where 
4ru,*= ee ee jarnrd BU ae) ed ee Aa 
ey! R 

corresponding to a surface distribution of sources, for the effect of A which 
corresponds to a surface distribution of doublets may then be deduced 
from it. it is also sufficient to assume that any line parallel to the v-axis 
cuts the body in at most two points. The integrand of (2.9) is exponen- 
tially small when & is large unless y—y’ and z—z’ are both vanishingly 
small so that if the line through (x, y, z) parallel to the x-axis does not 
meet the body at all u,* +0 as k-oo. If it does let one of the points 
of intersection be (€, y, z). Then in (2.9) B(y’, 2’) may be replaced by 
B(y, z) and taken outside the sign of integration and for the contribution 
to u,*, from the neighbourhood of (€, y, z) we must evaluate 


c y' dz 
pty. 2)| (e exp[k(x—2’—R)] . . . . (2.10) 


taken over the infinite plane lx’ +-my'-+-nz'=lé+-my--nz. 
Write 


(y y')(1 Pet lar cos O-+-nr sin @ Im(a—€)(1 [2)1/2, 
(z—z')(1—P?)1/2—Inr cos 6—mr sin 0+-In(a—€)(1— FP)” 
and (2.10) becomes 
rar 
eB ERG ae 


‘G d0 exp k{kl?(a—£)—(1—P)¥?r cos 8 
(le) 2) | 
I,(kr(1—P)"2) exp k[P(a—€) 
= (r+ Ae s)")*"]. (2.11) 


rar 
=rkB| [2+ P(@—é) zi fo 
—2rkB ifa>é ] 
=2rkB exp [2kl/?(a—€)] if sei 


(Doetsch 1943). Similar results hold for the second point of intersection 
between the line and the body and also for A. It is concluded therefore 


2) 


Z2A2Z 
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that (2.12) is the limiting value u* where now B is a function of y and z to 
be determined and é is that intersection of the line y=y’, 2=2' with the 
body nearer to (x, y, z). The results may also be written as 


* 
COE, as k +oo on the front side of the body, 
0x : 
(2.13) 
pac. as +0 as k-—oo on the rear side. 
0x J 
The boundary conditions on ¢ as k oo may now be determined. 
On the front side of the body 
0 0 
O= div ut=27 ge Se) +n [?k ut— SE oa ue i). . . (2.14) 
Oy oy az 


Since the first term is finite the last terms on the right hand must remain 
finite as k oo and so 


s~ 1 > OF ie nears (2.15) 
Then the normal component of u* — 0 as k + oo and therefore so also does 


the normal component of grad ¢. 
On the rear side when k is large, 


; =) (=) 
O=div u*¥=({ — aed ee Se ee 
(Fy ee Oz Jay ore 
Introducing a new coordinate B=2—E€ instead of x (2.16) implies that 
ss) (S) 

= = ==) a= 0 ee | 
(ay ae p (2.17) 

and since u*=-—grad ¢ on the body 


8%} OE 8 Pd OE Hd _ 


andy oy | oye conde koe Ge 
¢ being regarded as a function of x, y,z. Since ¢ is harmonic (2.18) 
reduces to 

ay Re ib a 0&0  - 0 (dd 


Thus on the rear side of the derivative ie along the normal is zero in 
the limit £ + a. 

In a two-dimensional flow (i.e. independent of z) (2.16) simplifies immedi- 
ately to v*=constant. Since 6=Uzx at an infinite distance downstream 
and u*=0 outside the wake (bounded by C) the integrated form of the 
equation of continuity shows that the constant must be zero. Hence 
0¢/0y—0 on the rear side of the body. 

If there is symmetry about the x-axis and 7 denotes distance from it 
a similar argument shows that d¢/dr=0 on the rear side. 

The boundary conditions which have been obtained above are all 
originally due to Oseen (1927). 
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At infinite Reynolds number therefore the ultimate flow is irrotational 
except in the wake which is bounded by the cylindrical surface C touching 
and extending downstream of the body and with generators parallel to 
the x-axis. This potential flow involves slip past the front side of the 
body which is reduced to rest by a boundary layer whose thickness is 0(v) 
from (2.12). This may be contrasted with the boundary layers in the full 
Navier Stokes equations whose thicknesses are 0(v1/2) in comparable cases. 
In the wake a flow with velocity components independent of « and which 
is rotational is added to the potential flow. This means that in the wake 
the fluid velocities do not tend to the uniform flow at infinity. 


§3. APPLICATION TO A SPHERE 

The results of §2 are now applied to the flow past a sphere, using a 
method originally given by Burgers (1921) for the flow past a circular 
cylinder. For that problem however Zeilon (1923) was later able to 
obtain the solution in closed form using a theorem of Hilbert. 

First we must find ¢. Let 7 denote distance from the axis of symmetry 
and write x=aR cos 0, r=aR sin 6,|@|<7 where a is the radius of the 
sphere. Then ¢ is harmonic and satisfies . 

eo¢/OR=0 if R=1 and w=|0|<n/2 3 1) 
dd/or =0 if R=1 and |6|<n/2 me 
and ¢~Uak cos 6 when Ff is large. 
Further we write 


= Ay , 


where the P,, are Legendre polynomials with »=cos@ and the A,, are 
constants to be determined in the course of the analysis. The third 
condition of (3.1) is automatically satisfied and so are the others if 


j=, — 2 A,,P,=0 —1xp<0 
n=0 
and ree (5-3) 
BTA. Pr 
M(u)=— sing Y — “0 Onl, 


n=0 n+ 1 
Following Burgers we solve these equations by the method of least squares 
minimizing 7 
x=| Lud | M(y)dp. eet 
=1 


For this purpose we make use of the following results which are either 
given or may be deduced from formulae in Modern Analysis (Whittaker 
and Watson 1927) 


‘0 1 
| P,P p= | fp \PP, dio 
—1 0 


if i—j is even and not zero ; and if 7 is odd j is even 
1 ue ] 2)P 1p ‘du= P,(0)P;'(0) : 
igen), ee ae alae Oe al 


0 ~ 
| Pee (3.5) 
=i 
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| 1 n+2 
— Of es le ee ee EE | ee eed | 
Then eae: ri 7 Ts) eats tat 
d 


(n—m)P,,(0)P,,'(9) 
(n—m)(n+m-+ 1)(m-+ 1)(n-+-m-+ 3) 


(eA ae 
HA ren 10+) : Ae eee ato 
ss gn (n—1)(n-+-2) +3: Pie ( ) 


For X to be a minimum 0X/dA,=0 for all A, and so we can set up an 
infinite set of simultaneous linear equations for the A,, with solution 


A,=—0-2645, A,=+0-4745, A,=—0-3201, Ap=+0-011 | 
A;=-10-100, Ag=—0-006, Ag=*40-054,0 Aj—4, 0-0042 eam) 
Ag=-+0-033, Az=—0-004, Ayp=—0-024 | 


while for higher values of n 


/ 


2 \1/2 2 ; 
A mo if n is even 


1/2\12 
(=) ; (=) if nis odd where ¢~10. 


Tf now a stream function Y=%)+7, is Introduced, defined by 
0f/or=ru, —dp%/Ox=rv 
where (uw v) are the components of u in the axial and radial directions, yy 


is the contribution from ¢ and %, is the contribution from u* it follows 
from (3.2) that 


5 Anbu lH) ay 
y= 2U a(1—p?)R?— Ua(1—p?) me n(n 1)R* +UaA,(1+p). (3.8) 
Further 4,—0 except inside C where 
7) 0 
a =— st on k=1 
; Op dy 
since u=0 there. Since ee =0 on R=1, 0<p<1 and J=0 on R=1, 
—1<p<0, J=0 on R=1, 0<p<l. 
Therefore since ys, is independent of a 
T i. AP 
b= Walp) + Vall) _Vadl+n,) . 8.9) 


where a?(1—p,”)=r?, if O<p,<1 and is zero otherwise. 

The streamlines are sketched in fig. 1 and their most striking feature 
is the slow inflow throughout the wake (r<a, «>0). As each streamline 
reaches C' its direction abruptly changes because of the discontinuity in 
u atr=a, «>0 and it moves outward to infinity. At large but not infinite 
Reynolds number there will be a shear layer along this surface and a 
boundary layer along the front side of the sphere. Further the wake 
will be closed and its length of order a2U/v so that the streamlines on which 
is is negative are also closed. 
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The pressure on the sphere is given by 


p=—p U0d/0x=pU? 2 A, P+ 4(u) fe EI), 
0 


Fig. 1 
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Streamlines of flow past a sphere at high Reynolds number according to Oseen’s 
approximation. The value of % on each streamline is marked on the figure. 


Fig. 2 


Pressure on the sphere. 


and is displayed in fig. 2. It is noted that the wake velocity at infinity 
which is in the direction of x, is given by p/pU in the range 0< | 6 |<7/2 
with r—asin@. The contribution D, from the pressure to the drag D on 
the sphere is 


af 4 ise 
D =—| 27a2 sin 6 cos OpdO=— = wpU*a* Ay. (3.11) 
p 4 3 
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so that the contribution from the pressure to the drag coefficient is 


In addition to the pressure a strict interpretation of Oseen’s theory 
means that there will also be a contribution from the skin friction, for the 
normal derivative of u* is 0(k) on the forward side and so when multiplied 
by the coefficient of viscosity a density of order one is obtained for the 
skin friction. With the full Navier Stokes equation the skin friction 
density would be 0(v/2) when v is small. Tomotika and Aoi (1950) have 
shown that on a sphere the skin friction contributes twice as much the 
pressure to the drag for all k and so in the limiting case k > oo we have 


D=3D,=1-06ra2U2p 
chat aeoie CHG: 


This result may be compared with the values of Cp, namely 2-85 at R=10 
and 2-08 at R=20, obtained by Goldstein (1929). 
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XXXII. 8Li Emission from Proton Induced Disintegrations at 950 Mev 


By B. A. Munir 
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[Received November 8, 1955] 


SUMMARY 
G5 plates exposed to 950 Mev protons were are ascanned and fifty *Li 


‘ Hammer stars’ discovered. They were divided into stars representing 
light and heavy nuclei disintegrations. In the case of light element stars 
the energy and angular distributions of the *Li fragment indicate that it 
is the residual nucleus after a light element breakup. For heavy nuclei 
stars the mechanisms of direct knock-on and evaporation can be excluded. 
The simplest hypothesis is that the fragments are ejected as part of the 
cascade process. 


§ L. INTRODUCTION 

A T-SHAPED track associated with a cosmic ray star was first reported by 
Occhialini and Powell (1947) and was thought to be due to a *Li fragment 
coming out of a nuclear disintegration and decaying at rest into two alpha 
particles and an electron. This was later confirmed by Franzinetti and 
Payne (1948) and Pickup and Voyvodic (1949). Alvarez (1950) pointed out 
that 3B should also give a similar ‘ Hammer track ’, though *Li emission 
should be more common due to a smaller potential barrier facing a ®Li 
nucleus. 

Studies of *Li fragments coming out of nuclear disintegrations have 
been made by several workers. Wright (1950) measured the cross 
sections for *Li production in light and heavy elements by bombarding 
them with 340 Mey protons and 190 Mev deuterons. He found a sharp 
discontinuity in the production cross sections in light and heavy elements. 
Hodgson (1951) analysed *Li fragments associated with cosmic ray stars 
and compared his results with predictions of the ‘ Evaporation theory ’. 
Titterton (1951) analysed the ejected *Li fragments from neutron induced 
stars at 150 Mev. He failed to see any °Li fragment associated with a 
heavy element breakup. ‘The results available so far, however, are not 
conclusive. In the previous analyses either the incident energy was not 
fixed or the direction of the incident primary was not uniquely defined. 
An attempt has therefore been made to investigate the ejected *Li frag- 
ments in proton induced stars at 950 Mev. The direction and energy of 
the primary being both fixed, it is hoped that the data may throw some 
further light on the mechanism of heavy fragment emission from nuclear 
explosions. 


+ Communicated by Professor P. B. Moon, F.BS. 
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§ 2. EXPERIMENTAL PROCEDURE 


G5 plates 400 microns thick were exposed to the 950 Mev proton beam of 
the Birmingham Proton Synchrotron, scattered out from the inflector tip 
(Moon, Riddiford and Symonds 1955). The beam was allowed to pass 
through the plates parallel to one edge and in the plane of the emulsion 
so that almost horizontal proton tracks were obtained. The plates were 
subsequently area scanned for stars under a x 10 objective and x 15 eve 
pieces. Each star with more than three prongs was closely scrutinized 
under a x45 objective and the prong number noted for stars having an 
incident minimum track in the general beam direction so that only proton 
induced stars were recorded. All black tracks from the proton induced 
stars were followed to their end in order to avoid losing very energetic 
SLi fragments. 

Stars near to the surfaces of the emulsion were not examined to mini- 
mize the loss factor due to the number of fragments passing out of the 
emulsion before coming to rest. 

A small percentage of the ‘hammer’ fragments may well be due to 
8B (Alvarez 1950) but in what follows it is assumed that all the “ hammer ’ 
fragments are $Li nuclei. 

The stars with ‘ hammer ’ tracks were divided into those due to heavy 
elements (Ag, Br) and those due to light elements (C, N, O). This was 
done because different mechanisms have been thought to be responsible 
for the emission of *Li from light and heavy nuclei (e.g. Wright 1950). 
It can be shown that stars due to light element breakup and having an 


associated ‘Li fragment can only have a prong number between four and 
seven. 


2C+1H SLi+3He-+ 21H (four prongs) 

160 4-TH sLi+-61H+3'n (seven prongs). 
In the case of seven-pronged stars, all the prongs except the ®8Li 
fragment must be due to singly charged particles. Thus all the stars 
with prong number greater than seven were classified as heavy element 
stars. In case of seven pronged stars, if any one prong other than the 
‘Li was identified as due toa doubly charged particle the star was included 
in the heavy element class. All stars with a short recoil were included 
in the heavy element category. The stars with prongs between four and 
seven, if not otherwise identified as being due to heavy nuclei were 
classified as light element disintegrations. Some uncertainty is possible 
in the above classification due to increase in prong number due to meson 
production but the energy distributions for SLi fragments from light and 


heavy elements (fig. 1(a) and 1(b)) suggest that the above classification is 
largely correct. 


§ 3. EXPERIMENTAL RESULTS 
A total of 25 200 stars of four or more prongs were examined and 50 
examples of the ejection of *Li fragments were found. Of these 28 came 
from heavy element disintegrations and 22 from light element disinte- 
grations. ‘The energy and angular distributions of the fragments are given 
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in figs. 1 and 2. In the latter figure the angle @ is the true angle in space 
between the ejected fragment and the direction of motion of the primary 
particle. The range-energy relation for the *Li fragments was calculated 
from that for protons as given by Rotblat (1950). 


Fig. | 
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Energy spectrum of the 8Li fragments ejected from heavy nuclei disintegrations. 
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Energy spectrum of the *Li fragments ejected from light nuclei disintegrations. 
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Fig. 2 
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(a) 
Angular distribution of the *Li fragments ejected from heavy nuclei 
disintegrations, with respect to the projected primary direction. 
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Angular distribution of the *Li fragments from light nuclei disintegrations. 
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The energy spectrum of the «-particles from the decay of the ®Be shows 
a peak near 3 Mev and agrees with the results of Baxter, Burcham and 
Paul (1950). All but six of the events had associated decay electrons ; 
the loss is probably due to observational inefficiency which arises from 
the finite thickness of the emulsions used. 


Heavy Element Disintegrations 


The energy distribution of the fragments (which has been corrected for 
loss in the usual way) shown in fig. 1 (a), has a cut-off at 8 Mev indicating 
that the potential barrier facing 8Li in a heavy nucleus is of this order of 
magnitude. The most energetic fragment observed was ejected with 
51 Mev. 

The angular distribution (fig. 2 (a)) shows a pronounced forward peak. 
the forward to backward ratio being 2-5: 1. 

Of the 25 200 stars examined, 21 800 correspond to disintegration of 
heavy nuclei (see Lock, March and McKeague 1955). Correcting for loss 
the number of ‘Li nuclei ejected from these stars is 30, corresponding to 
a production cross section of 1-1+-0-3 mb. 


Light Element Disintegrations 


The energy spectrum of the °Li fragments (fig. 1 (b)) has a peak between 
4 and 8 mev, the average energy is 12 Mev. The angular distribution 
(fig. 2 (b)) shows a forward to backward ratio of 2: 1. 

Correcting for loss the number of ®Li nuclei ejected from the 3400 light 
element disintegrations with four or more prongs is 22-5, corresponding to 
a production cross section of 0-5-- 0-2 mb. 


§ 4. INTERPRETATION 
8Li Emission from Light Elements 


The analysis of Titterton (1951) favoured the hypothesis that *Li is 
left as a residual nucleus in a light element breakup. Wright (1950) and 
Hodgson (1951) reached a similar conclusion. This recoil nucleus should, 
in general, come out with a small amount of energy and its angular 
distribution should be isotropic in the centre-of-mass system. 

The present results strengthen the above hypothesis. The fragments 
are emitted with a small amount of energy ; the angular distribution is 
anisotropic in the laboratory system, but at such high incident momentum 
the centre-of-mass velocity is not negligible. In the present experiment 
the incident momentum is 1650 Mev/c. In an encounter with a light 
nucleus, only a fraction of this momentum is transferred say, 300 Mev/c. 
This momentum given to an ‘SO nucleus corresponds to an energy of 
0-13 Mev per nucleon. This nucleus will be in a highly excited state after 
the impact and one can assume it to fly to pieces instantaneously, leaving 
a ®Li nucleus as a residue, which recoils with a small energy. For the 
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average recoil energy of 12 Mev, an isotropic angular distribution in the 
C-system will appear to be 1-8: 1 in the forward direction on consider- 
ations of the centre of mass velocity. This agrees with the observed 
distribution of 2:1 in the forward direction. Thus, the present results 
are compatible with the hypothesis that the *Li is a recoil nucleus from 
the light element disintegration. 


8Li Emission from Heavy Elements 

The simplest hypothesis to explain the emission of *Li fragments from 
heavy nuclei is that they are emitted during the evaporation process. 
Tf this is the case then their angular distribution should be isotropic, 
their energy spectrum of Maxwellian form and their frequency in agree- 
ment with that calculated by Le Couteur (1955). 

The experimental results do not seem to exhibit these characteristics in 
any striking way. The angular distribution is strongly peaked forward, 
with a forward to backward ratio of 2-5: 1 ; for heavy nuclei disinte- 
grations the centre of mass velocity is small and the correction to the 
angular distribution is negligible. The expected energy spectrum V(H)dE 
can be calculated using the relation 

(E—V) (HSV) 


Se meee 


Fe dE (Le Couteur 1950) 


N(E)dE= 


where the potential barrier V is taken to be 8 Mev and 7’ is the nuclear 
temperature. The mean number of tracks in heavy element stars with 
ejected ‘Li fragments is 9. On the average the number of shower and 
grey tracks per star is 2 (this agrees with the figures given by Lock et al. 
1955). It is reasonable to assume that 0-5 tracks per star are the low energy 
component of ‘cascade’ particles. Then the mean number of evaporation 
prongs per star is 6-5, including a SLi prong. This corresponds roughly to 
an excitation energy of 250 Mev and to an initial value of 7’~6 Mey. The 
energy spectrum calculated with these values is shown in fig. 2 (b) together 
with the spectra corresponding to 200 and 300 Mey excitation energies. 
It can be seen that the experimental distribution does not show good 
agreement with any of the calculated spectra. 

The frequencies of ejection of *Li fragments calculated by Le Couteur 
(1955) on the basis of evaporation theory are 0-0005 and 0-003 fragments 
per star for excitation energies of 150 and 250 Mev respectively. The 
experimental figure is 0-0014-++0-0002, in poor agreement with the cal- 
culated value for an excitation energy of 250 Mev. The agreement with 
the value for 200 Mev excitation energy is better but such an excitation 
energy is inconsistent with the average size of the disintegrations from 
which the fragments are ejected. 

Therefore, on the basis of the evidence discussed above, one may con- 
clude that it is unlikely that the ‘Li fragments are ejected during the 
evaporation process. 
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The results of Hodgson (1951) seemed to favour the evaporation mech- 
anism for *Li ejection from cosmic ray induced stars but his statistics 
were poor and he gave no quantitative measurements on the angular 
distribution. The results of workers on heavy fragments in general do 
not favour the evaporation hypothesis (see Perkins 1950, Camerini, Lock 
and Perkins 1952). 

The opposite extreme to the evaporation hypothesis is to assume that 
the *Li fragment is directly ejected by the incident proton. The maxi- 
mum energy of the *Li in such a collision would be 475 mev which is far 
im excess of the observed maximum energy of 51 Mev. This suggests 
that direct ejection does not often occur, if at all. It appears more 
reasonable to assume that the *Li fragment is ejected in some manner 
during the cascade process. Such a hypothesis is intermediate between 
the evaporation picture and direct ejection. The angular distribution 
is then readily explained since it is well known that the angular distri- 
bution of ‘ cascade ’ particles is peaked forward. The exact manner in 
which the *Li is ejected remains to be solved, as does the question of how 
such an aggregate of 3 protons and 5 neutrons keeps together even when 
ejected with energy up to 50 Mev. Presumably higher energy fragments 
are unstable and disintegrate immediately after leaving the nucleus. 


§ 5. Esection oF *Li From A DISINTEGRATING NUCLEUS 


During the scanning for *Li fragments, an event was discovered similar 
to a hammer but the alpha-particles were inclined at an angle of 140°. 
The event may be interpreted as the emission from a heavy nucleus of °Li 
which decays according to the scheme 


e€ 
9],i —_——+ *Be 


Due to simultaneous ejection of a neutron and two alpha particles from 
the excited Be nucleus, the alpha particles do not come out exactly in 
the opposite directions. The momenta of the alpha particles in this case 
are 122 mev/c and 88Mev/c and the angle between them 140°. The 
energy of the neutron can be easily calculated from this and is 3-5 Mev. 
Thus the disintegration energy of the ®*Be is 6-5 mev. Thus the event 
observed could be due to the decay of *Li to an excited state of °Be near 
8 mev, which subsequently gave two alpha particles and a neutron with a 
total energy of 6-5 Mev. : 

A similar event has been observed by Fry (1953) in plates exposed to a 
pe meson beam, which appeared to be due to a breakup of °*Be from 
the 6-8 mev level. 


n+ 24He. 


§ 6. CONCLUSIONS 


It seems probable that the explanation of Titterton (1951) is correct 
for ®Li emission from light elements. When a fast proton strikes a light 
nucleus, it transfers to it a part of its momentum, the remaining being 
carried away by fast knock-on particles or the incident proton itself. 


362 On the *Li Emission from Proton Induced Disintegrations 


The residual nucleus is in a highly excited state and almost spontaneously 
breaks up into pieces, not always charged, leaving “Li as a residual 
nucleus, which then recoils with a small amount of energy. The recoil 
energy should be in general higher, if a larger portion of the incident 
energy is transferred to the nucleus. Both the energy and angular 
distributions substantiate this assumption. The situation in the case 
of 8Li emission from heavy nuclei is not so clear. The * evaporation * 

and ‘ direct knock-on ’ hypotheses appear to be excluded by the experi- 
' mental data. The observations can best be explained by assuming that 
the ‘Li fragments are ejected during the ‘cascade’ process by some 
mechanism not yet understood. 
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ABSTRACT 


The single channel **Co gamma-ray spectrum and the spectrum in 
coincidence with K-x-rays have been studied with scintillation counters 
and a proportional counter and show no evidence for radiation other 
than the well-established 805 kev gamma-ray. Our results accord well 
with the decay scheme of Deutsch and Elliott. An upper limit of 
2-510 per disintegration was set on the intensity of 29 kev gamma 


radiation which, according to a recent suggestion, follows the decay 
of ®"Co. 


§ 1. IyrropuctTIoN 


THE intention of the experiments described in this paper is to resolve 
some of the uncertainties in the decay schemes of *8Co and *7Co which 
have arisen following a recent study of these isotopes by Cork, Brice 
and Schmid (1955) who find evidence for new gamma-rays in their decays. 

The situation regarding °8Co prior to their paper is embodied in the 
«lecay scheme presented in fig. 1 (a) which incorporates work reported 
by a number of authors. Livingood and. Seaborg (1941) detected 
positrons of maximum energy 500 kev (approximately) and gamma-rays 
of about 800 key. Deutsch and Elliott (1944) and Good, Peaslee and 
Deutsch (1946), using a beta-ray spectrometer, measured the end point 
of the positron spectrum and the converted electron groups from 507 and 
805 kev gamma-radiation. They also studied coincidences between 
Fe K-x-rays and gamma-rays and arrived at the branching ratio of 
positron emission and K-capture leading to the 805 kev level. Strauch 
(1950) determined the internal conversion coefficient for the 805 key 
radiation, as also did Cheng, Dick and Kurbatov (1954). Daniels e¢ al. 
(1952) measured the angular distribution of the 805 kev radiation from 
oriented ®8Co and on the basis of the decay scheme given in fig. | (a), 
were able to assign spins to the levels. 

The recent work of Cork et al. however favours the decay scheme given 
in fig. 1(6). They arrive at this conclusion from gamma-gamma 
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coincidence studies under conditions in which annihilation radiation alone 
is unlikely to give coincidences. Using scintillation counters they 
measured the pulse height distribution due to gamma radiation in 
coincidence with pulses attributed to the ‘ photo-peak’ of the 500 kev 
radiation and found that it was similar to the ungated spectrum. From 
this they have inferred the existence of a genuine 500 kev gamma-ray 


in the decay of ®*Co. 


Fig. 1 


O5 


O8l4 


FIG la FIG Ib 


(a) Decay scheme of Deutsch and Elliott (1944). 
(6) Decay scheme of Cork, Brice and Schmid (1955). 


For °7Co, in addition to the three well known gamma-rays of 14-4, 
123-0 and 137-4 kev (see Lemmer, Sagaert and Grace 1955 for references), 
Cork et al. have proposed further gamma-rays of 29, 99-8 and 700 kev. 


§ 2. SouRCE PREPARATION 


The cobalt 58 source was prepared through the ®8Ni(n, p)®*Co reaction 
by irradiating natural nickel in the fast neutron flux of the pile. In this 
way .effectively carrier free ®8Co can be separated without detectable 
amounts of other cobalt activities. This was then evaporated down on 
a sheet of thin cellophane. The cobalt 57 source was that used in earlier 
experiments and was deposited on 0-001in. aluminium foil. Its 
preparation has been described elsewhere (Lemmer et al. 1955). 


Decay of ®*8Co and 57Co 365. 


§ 3. CoBALT 58 


Two tests were made to distinguish between the two decay schemes. 
Firstly the ratio of the intensities of the 500 kev to 805 kev radiations 
was estimated from the pulse height spectrum in a scintillation counter : 
for decay scheme 1 (a) this should be 0-3 compared with 1-3. for 1 (6). 
Secondly the spectrum of gamma radiation in coincidence with x-radiation 
was determined : the x-rays arise almost solely from K-capture transitions, 
the internal conversion of the 500 and 805 kev radiations being known 
from previous evidence to be small. For decay scheme 1 (2) this should 
show only 805 kev radiation while for 1 (b) both 805 and 500 kev radiations 
should appear. 

The gamma-ray spectra have been taken with two Nal crystals both 
of 1jin. diameter and having lengths of 2 in. and } in. respectively. 
The spectrum taken in the 2 in. crystal is shown in fig. 2 (a). In order 
to determine the ‘ photo-peak ’ intensities of the 500 and 805 kev lines, 
the spectrum of a single gamma-ray of energy 845 kev from the decay of 
°4Mn was taken for comparison ; this spectrum is presented in the same 
figure after suitable scaling and adjustment in order to resemble more 
closely a gamma-ray of energy 805 kev. The ‘ photo-peak’ has been 
normalized to fit the ®°*Co ‘ photo-peak ’ and displaced by approximately 
5% in energy; the Compton spectrum has been decreased in height 
relative to the “ photo-peak ’ by 4%, this being the value predicted simply 
from the relative efficiencies of the two processes at these gamma-ray 
energies. The discrepancy, in the region below 400 kev, between the 
‘adjusted ’ 805 kev curve and the experimental points of figs. 2 (a), (b) 
can probably be attributed to the effect of external scattering; the 
54Mn source was further away from the detector than was the °5Co source. 
Even so it is possible to determine the ratio of the counting rates in the 
‘ photo-peaks ’ with an accuracy of 10%. 

To determine the true intensity ratio from this measured value, the 
total absorption of the two gamma-rays in the crystals was determined 
by numerical integration over the geometrical arrangement, the fractions 
of the counts in the photo-peaks being estimated from the shapes of the 
single line spectra ; a cross-check for this estimate is provided by the 
results of Maeder, Muller and Wintersteiger (1954). In this way the 
ratio of intensities of 500 kev to 805 kev radiation is found to be 0-29 
using the } in. crystal and 0-23 using the 2 in. crystal, giving a mean value 
of 0:26 with about 10°% accuracy. 

Coincidences were observed between Fe K-x-rays and gamma-rays, 
the latter being detected with the 2 in. crystal and the former, through 
a 0-0015 in. mica window, by a krypton filled proportional counter of 
the Maze type (West 1953). The source was placed between the two 
detectors, the sodium iodide crystal being covered by a thin lead sheet ; 
this prevented the 28 kev iodine K-x-rays or degraded gamma-radiation 
produced in the crystal from reaching the proportional counter. 

By means of coincidence and anticoincidence units coincidences were 
recorded only when the pulses in the x-ray channel lay in or near the 
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Fig. 2 
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(a) Gamma-ray spectrum from the decay of Co. Dots: experimental results. 
Full line: ‘ adjusted’ spectrum from the decay of °4Mn (see text). 

(b) Co gamma-ray spectrum taken in coincidence with the K-x-rays. 
Crosses : spectrum with background subtracted. Circles: background 
spectrum. Full line: ‘ adjusted’ °*Mn spectrum. 

(c) Spectrum of annihilation radiation from a ?2Na source. 

(d) K-x-ray spectrum from °°Co decay. The arrows show the bias limits for 
two * window positions ’. 
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peak. The coincidences were then made to open the kicksorter gate 
through which the pulses from the gamma-ray channel were fed. The 
resulting spectrum, after subtracting a calculated random background 
and a genuine background, is shown in fig. 2 (6) in which the genuine 
background spectrum is also plotted. Since the genuine background 
seems to account for all the coincidences in the 500 kev peak much 
importance was attached to its correct assessment. ‘ 

Figure 2 (d) gives the x-ray spectrum in the neighbourhood of the peak ; 
it will be seen that the peak is superimposed upon a slowly falling back- 
ground, thought to arise from detection of the positrons and gamma- 
radiation with incomplete absorption of the ionizing particle within the 
counter. Insertion of a thin aluminium foil between counter and source 
removed the x-ray peak completely and produced a reduction of the 
background in the neighbourhood of this peak, presumably by partial 
absorption of the positrons. The coincidence spectrum under these 
conditions was compared with that obtained with no absorber but with 
the x-ray channel ‘ window’ moved below the peak. The two spectra 
were not significantly different, since both positron-gamma and gamma-— 
gamma coincidences are expected to enhance the 500 kev peak relative 
to that at 805 kev. The statistical accuracy was not good enough to 
recognize small differences between these two sources of background. 
This comparison indicates that the coincidence spectrum due to back- 
ground in the proportional counter is unlikely to change appreciably 
when the “ window’ is moved just below the x-ray peak. The genuine 
background in fig. 2 (b) was obtained under the latter conditions. 

In fig. 2 (b) the adjusted Mn spectrum is also shown for comparison 
with the coincidence spectrum after removal of random events and 
background. It is seen that there is no significant peak due to 500 kev 
radiation and that certainly less than 10°, of the 500 kev radiation is 
in coincidence with x-radiation. 

These results cannot be reconciled with the decay scheme presented in 
fig. 1(b). Evenif the assessment of the background coincidence spectrum 
be completely in error, in which case up to 25% of the 500 key radiation 
could be in coincidence with x-radiation, only 6% of the 805 kev 
radiation could be part of the 500-805 kev cascade ; in fact we consider 
that our results give no significant indication of the presence of a cascade 
and further that not more than 24% of the 805 kev radiation could 
arise from such a cascade. 

In addition we have attempted to confirm more quantitatively the 
original decay scheme—fig. 1 (a). Our measurements on the intensity 
ratio of the 500 and 805 kev radiations are in agreement with this scheme, 
leading to a value of 13--1:3 °/, for the proportion of nuclei decaying by 
positron emission. 

As a confirmation of previous measurements on the positron spectrum 
we showed that coincidences occurred between the 805 kev and 500 kev 
radiations; sodium iodide crystals were used as detectors and one 
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channel was biased just below the 805 kev ‘ photo-peak’. It follows 
from the coincidence and single channel counting rates that the pre- 
dominant part of the 500 kev radiation is in coincidence with 805 kev 
radiation. 
§ 4. CoBALT 57 

We have examined the gamma-ray spectrum of °’Co with the krypton 
filled proportional counter in the energy range 9-48 kev. No evidence 
was found for any radiation other than that at 14-4 kev and in particular 
we can set an upper limit to the intensity of any 29 kev gamma-ray 
as 3%, of that of the 14-4 kev gamma-ray; this corresponds to an upper 
limit of 2-5 10-* per disintegration. 


§ 5. Discussion AND CONCLUSIONS 

The measurements on the decay of °*Co reported here are consistent 
with the generally accepted decay scheme (fig. 1 (a)) and not consistent 
with the decay scheme proposed by Cork, Brice and Schmid (1955). 
It is suggested that it may be possible to reconcile their data with the 
accepted decay scheme in the following way. They observed coinci- 
dences between two sodium iodide scintillation detectors one of which 
was gated on the 500 kev ‘ photo-peak’. The counters were arranged 
so that coincidences between annihilation radiation would not occur. 
Since they obtained in their coincidence spectrum a 500 kev peak having 
approximately the same intensity ratio to the 805 kev peak as occurred 
in the singles spectrum, they assumed that there must be a 500 key 
gamma-ray in cascade with the 805 kev gamma-ray. However, it is 
clear that there must be some pulses under the 500 key ‘ photo-peak ’ 
which are due to the Compton distribution of the 805 key radiation. 
Moreover the peak of the 805 key Compton distribution occurs at 
approximately 500 kev (see for example fig. 2 (a)). Therefore a 500 kev 
peak will appear in the coincidence spectrum due to coincidences with 
805 kev radiation, each 805 key quantum being associated with two 
annihilation quanta. It is easy to show that a result similar to that 
obtained by Cork, Brice and Schmid could be obtained by this mechanism. 

Our upper limit of 2-5 10-8 for the number of 29 kev gamma-rays per 
disintegration of °’Co makes it unlikely that the gamma-ray exists ; 
although Cork, Brice and Schmid give no value for its intensity it seems 
doubtful whether they could have observed a radiation as weak as this 
with a sodium iodide crystal. 
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XXxXV. CORRESPONDENCE 


Direct Observation of Static Nuclear Susceptibilities at Room Temperature 


By D. F. Evans 
Physical Chemistry Laboratory, Oxford 


[Received November 9, 1955] 


WHEN a substance containing nuclei with spin J is placed in a magnetic 
field, the nuclei occupy 2/-++1 energy levels, corresponding to the different 
orientations of the nuclear magnets with respect to the field. As a result 
of the Boltzmann distribution between these energy levels, they will not 
be occupied equally, and accordingly the substance has a resultant nuclear 
paramagnetic susceptibility y ) given by the expression 


we f+] 
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where N is the number of nuclei per cm?, and p is their magnetic moment 
(Smith 1953, Ramsey 1953). Nuclear susceptibilities at ordinary 
temperatures are extremely small ; thus for water at 20°c, since the proton 
has a spin of $ and a magnetic moment of 1-41 x 10-3 erg/gauss, the 
calculated nuclear susceptibility is 3-3 10-1! ¢.g.s. units, as compared 
with a diamagnetic susceptibility of 0-72x10-°. The only direct 
determination of a static nuclear susceptibility was made by Lasarew and 
Schubnikow (1930), who measured the susceptibility of solid hydrogen at 
4-22°, 2-18° and 1-76°K. At these temperatures, the nuclear susceptibility 
is an appreciable fraction of the normal diamagnetic susceptibility, and 
the observed (diamagnetic) susceptibility decreased as the temperature 
was lowered. | 

Nuclear magnetic resonance, which was discovered independently by 
Purcell et al. (1946) and Bloch e¢ al. (1946), involves the detection of 
transitions between these energy levels of a nucleus in a uniform magnetic 
field Hy, as the result of the application of a small radio-frequency field 
H, at the resonant frequency »H,/Jh. At resonance, the real and 
imaginary parts of the nuclear susceptibility y’ (dispersion) and  ,” 
(absorption) are both very much greater than the static susceptibility yo. 

The direct detection of the static susceptibility of nuclei at room 
temperature by the method of Lasarew and Schubnikow would clearly be 
extremely difficult. However, there is an alternative method available 
which involves the removal of this static susceptibility by means of a 
strong radio-frequency field at the resonant frequency. If the field is 
sufficiently strong, ‘ saturation ’ occurs, ie. as a result of the absorption 
of radio-frequency energy the populations of the energy levels are equalized, 
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and the nuclear susceptibility vanishes. The experimental problem is, 
therefore, to detect an increase of about 0-05° o In the diamagnetic suscepti- 
bility of a suitable substance when a relatively strong radio-frequency 
field at the appropriate resonant frequency is applied. 

The apparatus used was a modification of that previously constructed 
for work on photomagnetism (Evans 1955). Two thin glass tubes (~1 mm 
diameter, 9 mm long) containing solid hexaethylbenzene were suspended 
from a 20 quartz fibre in an inhomogeneous magnetic field (fig. 1). 


Fig. 1 
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BOTTOM VIEW 


Experimental arrangement. 


Hexaethylbenzene was chosen since it has a relatively high proportion of 
protons, and also in the hope that its relaxation time would not be unduly 
short (which would necessitate relatively high radio-frequency power 
levels to attain ‘ saturation’). The size of the silica balancing rod was 
such that, by suitable adjustment of the torsion head carrying the quartz 
fibre, the system took up a position similar to that shown in the figure, 
both with the field on and off. A Hartley oscillator supplied up to about 
100 volts across the radio-frequency coil. Since the samples are in an 
inhomogeneous field it is necessary to use a relatively broad Fedyen cy 
band to bring all the protons to resonance. This was ee by See 
the frequency (~40 times per second) over a range of 5-10°% by means of a 
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subsidiary condenser driven by a motor. At the resonant frequency the 
hexaethylbenzene becomes slightly more diamagnetic, and. accordingly 
the system swings in a clockwise direction (as seen from below). This 
motion was amplified and recorded as described previously (Evans 1955). 


Fig. 2 


> 
Increasing diamagnetism 


0 10 20 30 40 
Time (sec) 


The static nuclear susceptibility of the protons in solid hexaethylbenzene at 
room temperature. The arrow at the left indicates when the radio- 
frequency field was switched on, and that at the right when the field was 
switched off. The theoretical change in the diamagnetic susceptibility 
of the hexaethylbenzene (assuming complete saturation, which was 
probably not quite attained) is 0-049. The small zero drift is probably 
due to radio-frequency heating of the apparatus. 


A signal obtained in this way is shown in fig. 2. A mean frequency of 
about 30 Mc/s was employed, which was about that expected from a very 
rough measurement of the magnetic field. The effect of varying the mag- 
netic field and the frequency clearly showed that the signals obtained were, 
in fact, due to the nuclear susceptibility of the protons in the sample. 
Thus, when the magnetic field was increased, the effect virtually dis- 
appeared, but could be restored by increasing the frequency of the 
oscillator. Furthermore, if the oscillator motor was stopped, so that a 
very narrow frequency band was generated, only a negligible effect was 
observed. At frequencies and magnetic fields not corresponding to 
resonance, no effect was observed immediately after switching on the 
oscillator. However, if the oscillator was not quickly switched off, slow 
drifts were always obtained after a lag of about 10 seconds (probably due 
to radio-frequency heating of the system). As a result of the fairly long 
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period of the suspension, and the fact that the system was slightly over- 
damped, it is not possible to make an estimate of the relaxation time of 
hexaethylbenzene, although it must be less than about 3 seconds. 

It may be noted that the sensitivity of the above method of studying 
nuclear susceptibilities is not very much less than that of the conventional 
nuclear magnetic resonance. technique (only about 15 mg of hexaethyl- 
benzene were used, and the limit set by Brownian fluctuations had not 
been reached in: the present apparatus). However, it is unlikely to have 
any practical importance, since an inhomogeneous field is necessary and 
therefore the resolution is excessively low. Nevertheless, the experiment 
does demonstrate in a very direct manner both the existence of nuclear 
susceptibilities and the phenomenon of ‘saturation’. An even more 
direct method (suggested to the author by Dr. R. E. Richards) is to dis- 
pense with the radio-frequency field, and merely measure the magnetic 
susceptibility as a function of time immediately after switching on the 
main magnetic field. Some substances, such as lithium fluoride, have 
very long relaxation times, and accordingly the diamagnetic susceptibility 
would gradually decrease to a constant value as the nuclear susceptibility 
was excited by the magnetic field. An attempt to detect this effect was 
unsuccessful, due (at least-in part) to the instability of the apparatus 
immediately after switching on the field. A system with a carefully 
balanced taut suspension (as in a galvanometer) would almost certainly be 
more suitable. 


The author wishes to thank Dr. R. E. Richards for helpful discussions, 
and the University of Oxford for an I.C.I. Research Fellowship. 
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The Photoneutron Yields from ?°°Pb 


By J. D. Prentice and K. G. McNErILy 
Department of Natural Philosophy, The University, Glasgow 


[Received December 10, 1955] 


Tux relative photoneutron yields from ?°°Pb and natural lead have been 
measured at 22 Mev maximum bremsstrahlung x-ray energy, using the 
technique for determining (yn) yields developed by McNeill (1955). In 
the present work a specimen of lead enriched in the isotope of mass 208, 
isolated by Professor F. Soddy in his work on isotopes (Soddy 1916), and 
generously lent to us, was placed in the collimated beam of the Glasgow 
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22 mev synchrotron, and above a large Nal crystal. The photoneutrons 
created in the lead caused (n,y) reactions in the iodine of the crystal, 
producing #2°I, which decays by beta emission with a half life of 25 
minutes. The intensity of the activity thus produced was compared 
with the corresponding activity produced when the enriched lead was 
replaced by a geometrically identical specimen of ordinary lead. After 
subtraction of ‘ cosmic ’ background and the beam background produced 
when there is no target in position, the ratio of the activities produced by 
the enriched specimen and the ordinary lead sample is the ratio of the 
photoneutron yields in the two cases. This ratio was found to be 
Y enriched : Yordinary=1'09+0-16. The relatively high error is due to 
the small amount of enriched lead available, for the background was 
comparable with the activity produced by the photoneutrons from the 
lead samples. 

_ As the enriched sample contains 80° of ?°&Pb, with approximately 
10°, each of 2°*Pb and 2°7Pb, this result indicates that the ratio of the 
total photoneutron yield from lead 208 to the yield from ordinary lead is 
1-15+0-24. 

The opportunity has been taken to increase the accuracy of the order 
of magnitude measurement made by Reid and MeNeill (1953) on the 
22 mev yield for the photoproduction from natural lead of the 0-8 second 
isomer. It is concluded that the ratio of this yield to the 22 Mev yield of 
6&Cu produced by ®Cu(yn)®Cu, written Yoo/ Yay, 18 0-27, which is 
somewhat higher than the original estimate. 

In conjunction with other known photoneutron yields, the present 
results can be used to find the proportion of photoneutron events in 2°°Pb 
which lead to the isomeric state of 2°’Pb, that is, Yoo7+/Yoo,- For, as was 
shown by the use of enriched isotopes (Reid and McNeill 1954) the 0-8 
second isomer is produced by a photoneutron reaction in ?°8Pb, and the 
22 mev photoneutron yields from ®Cu and ordinary lead are weil 
established. 

Then (¥ so7«/ ¥ 297) (¥ soze/ Y sou) X (Yercu/ Y pp) X (Vpp/Y 207) where the 
subscripts refer to the product of the photoneutron reaction, and Y p, refers 
to the neutron yield from natural lead. 

As stated above, Y97«/Yuoy has been found to be 0-27, and Yug,/Y pp, 
is 2-8/27 (Katz and Cameron 1951, and Montalbetti et al. 1953). The 
present work shows that the final term is 0-87-+0-18, giving a final 
result of Yogo7«/¥o9,=0-025. It may be concluded that about 25 in 
every 1000 photoneutron events induced in ?°Pb by bremsstrahlung 


X-rays of 22 Mev maximum energy lead to the 0-8 second isomeric state 
iy), 207P 
hier dl ta a 
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XXXVI. REVIEWS OF BOOKS 


Meteors. Edited by T. R. Karsmr. (London: Perg Press 
cn ae ree ( ergamon Press Ltd.) 


Tuis work consists of 39 papers by authors from some ten countries, presented 
at a symposium on meteor physics held at the Jodrell Bank Experimental 
Station in July 1954. It forms also a special supplement (Vol. II) to the 
Journal of Atmospheric and Terrestrial Physics. 


Static and Dynamic Electron Optics. By P. A. Srurrock. (Cambridge : 
University Press.) [Pp. x+240.] Price 30s. 


A BRIEF explanation of the title is needed. Sturrock deals with mathematical 
techniques which originated from the study of electron microscopes where 
electrons move in static fields. But the scope of these techniques goes far 
beyond electron microscopy. In particular, in the second part of his book, 
entitled Dynamic electron optics Sturrock extends them to accelerators, with 
time-varying orbits: he succeeds in bringing two subjects under one roof. 

Indeed, the book teaches the mathematical craft of tracing charged particles 
through quite general electromagnetic fields in a meaningful manner, i.e. giving 
the relevant information in each case, such as focus, quality of picture, stability, 
beam width, etc. Thus the book is useful to anyone working, say, in electron 
tube design, electron microscopy, mass-spectroscopy, high-energy accelerator 
physics, cosmic-ray physics. It ought to be consulted before calculations are 
begun on the instrument or situation in hand. 

For Sturrock’s methods are analytical, general and powerful. He exploits 
fully the beautiful theories of Hamilton, Lagrange and Poincaré and takes 
these beyond the realm of aesthetic mathematical pastimes and_ tripos 
exercises. In many respects, the book is an apt successor to Synge’s Geo- 
metrical Optics and De Broglie Waves of the same series of monographs. 
The variational method, ‘ perturbation characteristic functions * and invariants 
are the basic tools. In all cases the establishment of an axial (unperturbed) 
ray is followed by the discussion of the paraxial (linear) approximation which 
is then improved upon by the calculation of aberrations or * non-linear ’ effects. 

In this age of mechanization, it is gratifying to find that Sturrock has banished 
numerical drudgery from his calculations. By demonstrating the power of the 
analytical method, he has rescued the subject from being devoured wholesale 
by computing machines to whom orbit integrations are very digestible fare. 
Indeed, a problem would have to be of a very special and highly non-linear 
kind before Sturrock’s perturbation methods failed and had to be complemented 
by computing effort. The successful handling of non-linear effects in accelerators 
in the last chapter of this book demonstrates the superiority of his methods. 

He is able to generalize static electron-optical methods to the ‘ dynamic ’ 
accelerator problem with such ease because, throughout, he observes relativistic 
accuracy—a comforting thought in all applications—so that time, to him, is 
‘just another coordinate”. In fact, the only remaining distinction between 
‘static’ and ‘dynamic ’ electron optics lies in the different interpretations of 
the term ‘ focusing’. The book is not ‘easy’. The arguments are condensed 
and the mathematics, although not beyond the range of an average physicist, 
is not ‘elementary’. It therefore remains a specialists’ book. One has the 
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feeling that a little expansion here and there would render the book readable 
to a wider public. Some changes in notation are called for : the letter * p 
is used in too many senses, albeit in different type. Some of the formulae 
could be symmetrized : equations which are in principle covariant (especially 
the fundamental equations) are not manifestly so. 

But these are minor criticisms which should not deter the serious student of 
the subject from reading this highly commendable book. O. B. 


Thermodynamics and Physics of Matter. (Volume I of a Series on High Speed 
Aerodynamics and Jet Propulsion.) Edited by F. D. Rosstnr. (Oxford : 
University Press; Princeton: University Press.) [Pp. xviii+812.] Price 
100s. ($15.00). 

Hiau speed flight and propulsion have now advanced to a stage where the 

designer or research worker needs to understand and apply a very wide range 

of information, embracing applied mechanics, physics and chemistry. This 
situation has prompted an eminent American committee to begin the com- 
pilation of a comprehensive series of twelve volumes designed to present all 
the fundamental information conceivably needed by modern flight and pro- 
pulsion engineers and, in the later volumes, to indicate its utility and appli- 
cation. This encyclopaedic aim has necessitated the division of not only the 
series but also each individual volume among serveral specialist authors. One 
danger of this arrangement is that the marshalling of the sections in logical 
order and with uniform notation without omissions or undue repetition may 
prove too difficult. On the whole the material has been well organized in the 
volume under review. Another danger is that the specialist authors may lack 
the space or the enthusiasm to treat their topic fully and instead refer the 
reader to publications elsewhere. When this happens, as it does occasionally 
in Volume I of the series; one begins to doubt the desirability of such com- 
pendious publications, since they are evidently not a substitute for specialist 
works. If space is at a premium in these volumes, the allotment of space to 
introductory elementary thermodynamics can only be attributed to the need 
to counteract the wide distribution of engineering thermodynamical texts 
which treat the subject in a form ill-suited to more advanced development. In 
contrast the sudden, if inevitable, introduction of statistical mechanics with- 
out preamble in the first section on thermodynamic fundamentals is not 
calculated to encourage the reader with the normal engineering background. 

A more thorough study of the implications of statistical mechanics occurs in 

later parts of the book. 

There are comprehensive chapters on gaseous, liquid and solid state theory 
and on transfer processes and reaction rates in the fluid states, and a very 
detailed study of the behaviour of fluids near the critical point. It is obvious 
by this stage that the book is aimed at a wider public than workers on high 
speed flight and propulsion. The section by Beattie on the properties of real 
substances and mixtures will however prove invaluable to aerodynamicists 
and engineers. The relaxation phenomena important in severe shockwaves 
and elsewhere receive very adequate treatment. A section on gases of such 
low density that continuous fluid theory is inapplicable paves the way for 
aerodynamics at high altitude, presumably to be discussed in later volumes. 
The volume ends with a brief study of the thermodynamics of irreversible 
processes. Nearly all the chapters end abruptly and would have been im- 
proved by the addition of general remarks summarizing the important points 
and indicating their relevance to high speed flight and propulsion. J. A.S 
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Spectroscopy at Radio and Microwave Frequencies. By D. J. E. Ineram. 
(London : Butterworths Scientific Publications.) [Pp. xii+332.] Price 45s. 


ALTHOUGH spectroscopy at radio and microwave frequencies was almost 
unknown ten years ago, the subject has developed rapidly and this review of 
the main applications and techniques is welcome both to the research worker 
in the field and to the non-specialist who requires an idea of methods and 
possibilities. 

After discussing the relation between microwave and ordinary spectroscopy 
and indicating the scope of the subject, the author devotes the first three 
chapters to the production and measurement of microwave power, frequency 
measurement, waveguide techniques, and the design and construction of 
different types of microwave spectroscope. 

The following four chapters are concerned with the theory and the experi- 
mental results for the three main branches of the subject, gaseous microwave 
spectroscopy, paramagnetic resonance and radiofrequency spectroscopy. In 
each case the theory is stated in a clear and concise manner without detailed 
mathematical treatment and well illustrated by reference to the experimental 
results which are given in considerable detail. 

The final chapter is devoted to a very useful discussion of the applications 
of microwave and radiofrequency spectroscopic techniques. The inclusion 
of these topics and the large number of references to original papers at the 
end of each chapter should make this book of considerable interest to those 
engaged in both pure and applied research. KAS: 


[The Editors do not hold themselves responsible for the views 
eapressed by thear correspondents. | 


ADDENDUM 


Two omissions in the bibliography of ““ The Thermal Instability of a Fluid 
Sphere Heated Within”, Phil. Mag., 46, 1310, have come to the attention 
of its author. They are: 


Maturkuar, 8. L., 1937, Proc. Ind. Acad. Sci., 7,,.No. 1, Sec. A. 
JEFFREYS, H., and Jerrreys, B. 8., 1950, Methods of Mathematical Physics 
(Cambridge : University Press), p. 442. 


ERRATUM 


The Magnetic Susceptibilities of some Diamagnetic Alloys, by W.G. Henry 
and J. L. Rogsrs, 1956, Phil. Mag., 1, 246. 


The denominator in eqn. (4.9) should read 4:mc? not 8mec?. 


